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ABSTRACT: newly constructed single-opening marina was evaluated determine 
flushing characteristics comparison hydraulic model and mathematical model 
with field studies dye releases. Both models simulated well-mixed conditions and 
agreed well with each cther. The dye study was conducted mixing rhodamine-wt 
throughout the marina flood tide and monitoring dye concentration versus time 
over several tidal cycles following release. Over the initial 7-10 hours the field 
studies, ail three methods agreed well; field results thereafter showed three six times 
less dye than model predictions. For marinas simple geometry with single openings 
and width-to-length ratios similar the marina studied here, flushing estimates may 
conservatively approximated elementary mathematical methods. 
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ABSTRACT: Empirical distributions wave heights derived from field data differ 
from the theoretical distribution wave heights given the Rayleigh law. 
principal cause this discrepancy suggested that the conventional theory does 
not any way account for wave breaking. Based the assumption that the sea 
surface characterized with narrow-band spectrum, the probability distribution 
wave heights limited breaking, due Stokes-type limiting steepness deep 
finite water depths, derived. The breaking-limited distribution provides probability 
discription which compares with two well documented field observations more 
favorably than the Rayleigh law. 


REFERENCE: Tayfun, Aziz, “Breaking-Limited Wave Heights,” Journal the 
Waterway, Port, Coastal and Ocean Division, ASCE, Vol. 107, No. WW2, Proc. Paper 
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ABSTRACT: test pipe rig was designed, fabricated, and then placed the ocean, 
offshore from Honolulu. The pipe diameter was inches, its length 17.5 feet, and the 
water depth was feet. During four seasons swell from Southern Hemisphere 
winter storms, data were obtained for the forces exerted 40-inch-long portion 
the test pipe. Four separate angles between the swell fronts and the test pipe were 
involved: 0°, 15°, 35°, and 52°. the measured forces were combined with concurrently 
measured wave-induced kinematic values derive horizontal and vertical pipe force 
coefficients use offshore pipeline designers. 


REFERENCE: Grace, Robert A., and Zee, Gabriel T.Y., Forces Rigid Pipes 
Using Ocean Test Data,” Journal the Waterway, Port, Coastal and Ocean Division, 
ASCE, Vol. 107, No. WW2, Proc. Paper 16272, May, 1981, pp. 71-92 
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ABSTRACT: The various steps involved the estimation design wave 
are reviewed. These include the following: methods data collection, the use 
plotting formula, the selection suitable distribution and its parameters, the plotting 
confidence bands about the best-fit line, and the selection design wave 
corresponding prescribed return period encounter probability. The probability 
distributions commonly used, including the log-normal and Extremal Types II, and 
III distributions (Gumbel, Fretchet, and Weibull respectively), are described, and their 
properties are summarized tabular form for ready reference. 


REFERENCE: Isaacson, Michael St. Q., and MacKenzie, Neil G., “Long-Term 
Distributions Ocean Waves: Review,” Journal the Waterway, Port, Coastal and 
Ocean Division, ASCE, Vol. 107, No. WW2, Proc. Paper 16277, May, 1981, pp. 93- 
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FLUSHING STUDY SOUTH BEACH 
MARINA, OREGON 


Richard 


INTRODUCTION 


Recent increases recreational and small commercial craft activities have 
resulted the construction many new marinas. Local, state, and Federal 
governments must evaluate applications for marina construction permits. Little 
information exists ecological impacts marinas construction events 
such dredging and spoil disposal. 

This report concerns one aspect the marina permit evaluation process: 
water quality impacts related marina circulation and flushing efficiency. These 
physical properties vary with the tide range, water density, and physical 
dimensions marina. Water quality is-affected the degree flushing, 
and sediment redistribution currents. Detrimental water quality can determine, 
e.g., the fate migrating juvenile fish and benthic organisms (5,7). 

Ecological studies marinas are few. The most comprehensive have been 
performed Marina del Rey, Calif. (3,15,16,17,22). Slotta and Noble (21) 
analyzed the use benthic sediments indicators marina flushing several 
Pacific Northwest marinas. Puget Sound marina water quality studies have also 
been conducted (10, 11, and unpublished work Yearsley). 


Marina 


Most mathematical studies marina circulation and flushing have been 
concerned with vertically well-mixed waters one and two dimensions (1,2,4,6, 
18,20). 

Hydraulic model studies small harbors provide alternative method 
assessing flushing ability, although not without deficiencies related scale 
distortion. Several hydraulic model studies Pacific Northwest marinas have 
been conducted Nece and Richey and their associates the University 
Washington (9,12,13), and Slotta and others Oregon State University 
(1,20,21). 

Environmental Protection Agency, Marine Div., Environmental Re- 
search Lab., 200 S.W. St., Corvallis, Oreg. 97330. 

Note.—Discussion open until October 1981. extend the closing date one month, 
written request must filed with the Manager Technical and Professional Publications, 
ASCE. Manuscript was submitted for review for possible publication July 16, 1980. 
This paper part the Journal the Waterway, Port, Coastal and Ocean Division, 


Proceedings the American Society Civil Engineers, Vol. 107, No. WW2, 
May, 1981. 
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Because the small size most Pacific Northwest marinas, use numerical 
models employing finite difference grids finite elements not always practical. 
order utilize the numerical method its best advantage, very small grids 
would need used; for finite difference analogs employing explicit solutions, 
this would turn require small time steps less than the grid size, Ax, divided 
the speed shallow water wave, gh, for one-dimensional simulation. 
For 328-ft (100-m) grid size in, e.g., 20-ft (6-m) deep marina [depth 
3.3-ft tide amplitude], the time step would need less than sec. 
Smaller grid sizes greater depth, require proportionately smaller time steps 
and increased computer time. 

The Oregon marinas examined Oregon State University used Froude scale 
models, did those constructed the University Washington. All the 
model studies (except that reported Ref. 13) were site-specific. 

One advantage hydraulic model over numerical one that small-scale 
operations relating mixing can readily observed and photographically 
recorded. These phenomena approximate mixing processes which take place 
the prototype. attempt has been made study vertical exchange processes 
either the numerical hydraulic models. Rather, time and space averages 
were taken; the water column assumed well-mixed the vertical. 
This assumption prevents the reproduction vertical horizontal convection 
currents. 


Surveys 


Dye releases have been used determine flushing rates Florida finger 
canals (2,8,23). Slotta and Tang (20) released dye Oregon’s Chetco estuary 
boat basin, and compared results with hydraulic and finite element model. 
Discrepancies between field and hydraulic model results were due the difficulty 
obtaining the proper density differences between dye and the receiving water 
the hydraulic model. However, depth-averaged concentration versus time 
curves were similar. 

the field experiments described, dye tracer was distributed 
throughout the marina during the middle stage flood tide. The last two 
hours the flood tide were used permit the dye continue mixing. 
initial average dye concentration, C,, was achieved maximum high tide; ideally 
this concentration remains constant the following ebb tide, while dye mass 
decreases. Assuming return dye the next flood tide, the mass 
dye the marina remains constant while the concentration decreases with 
increase volume during the flood. 

These assumptions can expressed follows, for dye mass, concentra- 
tion, with increasing volumes represented decreasing volumes 
constant volumes for flow, and volume, 

Note that assumed that the dye uniformly mixed throughout the basin 
the time initial high tide. 
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For the case direct fresh-water inflow, the volume water present 
marina can evaluated terms mean tide level volume, V,, and 
variations about the mean. For mean low tide volume A,R/2, 
which mean marina surface area; and the tidal range from mean 
iow mean high water. 

water ebb tide carrying with mass, (for the first ebb), concentration, 
the end the ebb the mass is, the preceding 


M, = M,- C,(A,R) 


the following flood tide, the concentration decreases; the end the 
flood 


which depth mean low water level, while the mass unchanged. 
There stepwise decrease both and however the time change between 
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FIG. 1.—Yaquina Estuary Entrance and South Beach Marina Place Map) 


constant values gradual rather than abrupt, and and are out phase 
90°. follows that the concentration after the ith flood cycle 
This relationship can explored through analysis first-order 
differential equation. The rate change concentration assumed proportional 
the concentration present: 


the frequency tidal period, and phase angle. 
Solving numerical differentiation 


= 
WASH. 
aff 7 
= 
a a 
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where the subscripts time; and the time step. The are any given 
time units, not necessarily tidal cycles, with the restriction that ebb cycle, 
C,. The equation can also solved direct integration. 

South Beach Marina entrance South Beach marina 
1.5 nautical miles (2.8 km) upstream the end the north jetty Newport, 


FEET 


FIG. 2.—Exchange Coefficient South Beach Marina Hydraulic Model (19); 


Oreg. (see Fig. 1). Approximate marina dimensions are: length, 1,574 
(480 m); width, 623 (190 m); depth, mean tide level 13.4 


(4.1 m), resulting mean tide level (MTL) volume, 1.34 
10° The mean tide range, 6.0 (1.83 m), resulting mean 
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The entrance width 157 (48 m), resulting mean cross-sectional entrance 

The marina designed hold 600 boats. the time the field work, 
boat slips piles were place. Dredge spoil was disposed the east 
bank the marina. The marina breakwater jetty restricts rapid flowthrough. 

South Beach Marina Hydraulic Model hydraulic model study 
South Beach marina was made Richey and Skjelbreia (19). model 
used 1:10 distortion ratio; Froude scaling was employed 
with the following scale ratios: horizontal length 1:480; vertical length 
1:48; velocity 1:6.93; time 1:69.3. The prototype tidal cycle was taken 
12.4 which required 10.74 min reproduce the model. Sinusoidal 
tides were used for 0.9-m, 1.8-m, and 2.7-m ranges. Model water density was 
boat slips and pilings and wind stress were not modeled. 

20% solution rhodamine-WT diluted 1:100 was used tracer. Thirty 
milliliters dye were completely mixed the model high water location. 


> 
= 


FIG. 3.—Relative Dye Concentration Versus Tidal Cycle, South Beach Marina 
Hydraulic Model (19) 


Turner Model 110 fluorometer was used determine relative concentration 
the end four tide cycles. 

This procedure was augmented time-lapse photography the loss 
water-soluble dye with time. The photos were analyzed densitometer 
obtain relative concentration versus time. Exchange coefficients, defined 
were determined from both the dye and densitometer 
data. Here, concentration the ith tidal cycle; and initial concentration. 
Exchange coefficients thus calculated are shown for the 5.9-ft (1.8-m) tide 
Figs. and Fig. shows values the end four cycles. These are 
averaged values and range from about 0.15 the southwest corner 0.5 
the northeast corner. Low values correspond relatively low flushing, while 
higher values indicate greater flushing. This shown Fig. where individual 
C,/C, values are 0.35 for station and 0.1 for station This corresponds 
0.23 and 0.44 for stations and respectively, which the 


Scale Sampling Stations 
x4 
0.8 
0.6 
0.2 
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range the values given the previous figure. 

The main conclusions the hydraulic model study were follows: (1) The 
basin had good hydraulic characteristics except for poor exchange the south 
corners; (2) good exchange due strong currents the main channel passing 
the marina entrances, which can introduce turbulent eddies flood tide and 
prevent recirculation the ebb; and (3) the currents past the entrance improve 
the exchange coefficient about 25% over that where the entrance situated 
quiet bay. 

Field Float Studies.—Weighted poles 6-ft, 8-ft, 10-ft, and 12-ft (1.83-m, 
2.44-m, 3.05-m, and 3.66-m) submerged length were released and followed 


ww2 


Start 


and Finish 
Dye Re- 
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200 
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FIG. 4.—South Beach Marina Field Sampling Stations, 1978-1979; Dashed Line 
Approximate Dye Release Track 


January 17-18, 1979, conjunction with dye release; pole positions were 
fixed sextant. Easterly winds ranged from knots (1.5 
m/s) the first survey day which was conducted ebb tide; 10-knot 
winds determined the main direction all pole trajectories during this ebb 
cycle although northwest water current component toward the entrance was 
present. Tide decreased from maximum high water (LHW)] height 
6.9 ft-3.9 (2.1 m-1.2 m), about midtide. Maximum pole velocities observed 


for the 6-ft and 8-ft poles were 0.8 fps and 0.9 fps (26 cm/s and cm/s), 
respectively. 


> 


ww2 SOUTH BEACH MARINA 


Winds were less than knots from the west the January survey. Tides 
increased from 6.2 ft-7.5 (1.9 m-2.3 [higher high water Although 
this study was only conducted for short time, pole trajectories show that 
surface water near the entrance had net outward direction, while the deeper 
layer showed inward motion. Thus, there was some flushing the incoming 
tide the upper layers. 

Field Dye was released for about over the entire 
marina (see Fig. starting about before slack water flood tide. Visual 
observations from the United States Highway 101 bridge Newport and 
light plane did not reveal any obvious high low surface dye patch concentrations. 

Fifty-one (23 kg) 20% dye was diluted with gal (190 seawater 
approximate receiving water density from 55-gal (210-L) drum; was 
discharged about 0.8 gal/min L/min). Two methods release were used. 
the first survey, 15-ft (4.6-m) hollow aluminum pole with horizontal 
discharge tube the bottom was raised and lowered the dye was released. 
The outboard motor vessel made right angle paths throughout the marina. 
the second survey, garden hose was towed 45° from the surface 
bottom behind the boat throughout the water column; discharge ports 
the hose about 1.5-ft increments allowed for more even vertical distribution 
the dye than occurred the first survey. 

the locations shown Fig. hourly samples the surface, middepth, 
and about 1.5 from the bottom were drawn from continuous flow hose-pump 
arrangement. Samples were taken hourly intervals for the first and 
mid- and maximum high and low tides thereafter until background levels 
were approached. 


Samples were analyzed board Turner 111 flow-through fluorometer 
fitted with 546- excitation and 590- emission filters. Samples for analysis 
the laboratory were collected 125-ml screw-cap containers after local 
equilibrium was reached indicated the field fluorometer. The Turner Designs 
fluorometer was used the laboratory. Frequent field calibration was made; 
laboratory standards were used before and after each run, values reported here 
are from the laboratory analyses. 


The 1978 study was similar that 1979 except for the following differences: 
(1) Average tide ranges were 7.9 and 4.9 (2.4 and 1.5 m), respectively; 
(2) water column density was different, with greater rate change density 
with depth for the 1978 study; and (3) the method introducing the dye was 
different. 

Stratification was slight each study, the main differences being due 
temperature. Salinity and temperature-depth profiles showed gradual decrease 
temperature from and increase salinity from 31.6 
32.1 from surface bottom for 1978. 

1978 Survey.—Attempts distribute dye evenly throughout the water column 
were not completely satisfactory. Fig. shows surface, mid, and bottom dye 
concentrations station (refer Fig. 4). The surface and middle concentrations 
were nearly equal about 1400 September 15, 1978, but the bottom sample 
was initially quite low; all samples approached equality after about 0600 


= ea. 
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September 16, 1978. Surface concentrations for all six stations over the duration 
the experiment showed similar concentrations except for station which 
decreased rapidly about 1800 the first day. Bottom concentration with 
time showed station having lower concentrations than all the others 
2200 (September 15) when, again, the concentration lines merged. Stations 
and showed several spikes which were not apparent the other values. 


CONCENTRATION (ppb) 


00 (PST) 
— 9/15/78 9/16/78 


FIG. 5.—Surface, and Bottom Rhodamine-wt Concentration, parts per 
billion, South Beach Marina, September 15-16, 1978 


@ 


(2) 


Hours After High Water Slack 


FIG. 6.—Relative Dye Concentration (C,/C,) Averaged Over Depth and Model 
Computations Versus Time, South Beach Marina, September 15-16, 1978 


Fig. shows values averaged for all stations. The modified tidal 
prism and numerical model results are also shown. Considering the variability 
the dye patches, the first are reasonable agreement. rapid decrease 
concentration shown from hours 7-12 the models and the field data. 


x 
' 
BOTTOM: 
/ 
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The latter, however, continue show decreases about hour 18. 

The 1978 field data can also compared with the hydraulic model results 
for the 8.8-ft (2.7-m) range tests. Exchange coefficients the hydraulic model 
ranged from Interpolation between the 6-ft (1.83-m) and 9-ft (2.74-m) 
give 0.44 for the 8-ft (2.4-m) range found the field. The exchange 
coefficient based the prism method 0.5 for the 8-ft range. This gives 
values 0.55, 0.30, ... alternate high tides starting with C,/C, 
Field, hydraulic, and mathematical model results all clearly show rapid flushing 
for this tide range. The hydraulic and mathematical model results are essentially 
equivalent, but underestimate the dye removed when compared with the field 
study. 

1979 Survey.—For the 1979 survey, the sampling station location was essentially 


the same before except that station the 1978 survey was eliminated 
(Fig. 4). 


Ci/Co 


Hydraulic Mode! Averages 
(See 


Over Depth 


oO 

— 


Hours After High Water Slack 


FIG. Dye Concentration Averaged Over Depth and Model 
Computations Versus Time, South Beach Marina, January 17-18, 1979 


Plots (not shown) all stations the surface, middle, and bottom show 
better initial mixing than 1978. Station the northeast corner, initially 
showed higher concentrations than the other stations but quickly merged with 
the rest. All data showed increase C,/C, ratios from about hour 3-7 
(ebb cycle), while values station indicated recirculation south along the 
west side the marina. 

The hydraulic model results (19) show average interpolated exchange 
coefficient 0.3 for the 6-ft range. The high tide C,/C, values station 
the model and field are close; there initial rapid decrease concentration 
about C,/C, 0.2, followed gradual decrease 0.1 the fourth 
tide cycle. Fig. shows the average C,/C, results for the hydraulic model 
taken from Fig. and the mathematical model values. 


1.2 
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curves merged after the second flood tide. 


Comparisons Figs. and show that normalized concentration-time curves 


were quite similar with respect the sharp decline during the first flood cycle. 


The 1979 curve shows steeper exponential decrease toward background 
concentration. 


Summary 


South Beach marina has single entrance and uncomplicated geometrically; 
has free exchange with the main navigation channel where rather large currents 
develop during ebb and flood tides. 

Mathematical and hydraulic models results agreed well when exchange coeffi- 
cients were averaged over the entire hydraulic model. 

Flushing efficiency near the marina entrance about twice great the 
inner harbor indicated the hydraulic model over four high tide cycles. 
Field results show similar spreads although variations are more extreme among 
stations. 

Neither the hydraulic nor mathematical models successfully reproduced early 
flushing events; both underestimated the first flood decrease concentration 
about 30%-40%. Thereafter, the predicted and observed curves parallelled 
each other, although the model-predicted curves remained higher. terms 
pollutant concentrations, both model results were more conservative. 

Hydraulic model studies successfully predicted that the South Beach marina 
would have satisfactory exchange because the strong currents moving past 
the entrance. These entrance currents undoubtedly provide significant transfer 
processes through vortex motion and gyre generation. unlikely that marina 
similar dimensions would well flushed were sited less active 
environment. 

The mathematical and hydraulic model comparison raises doubt the 
need for hydraulic model studies small marinas one can content with 
conservative predictions. For marinas similar dimensions similar locations, 
good approximations flushing efficiency can made matter minutes 
using simple box-model assumptions. However, these results cannot extrapo- 


lated other marinas with multiple entrances, significantly different width-to- 
length ratios, environmental settings. 
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The following symbols were used this paper: 


cross-sectional planform area; 
concentration; 

mean tidal depth; 

exchange coefficient; 

gravity; 

maximum depth, shallow water wave; 
tide cycle; 

marina length; 

dye mass; 

tide range; 

tidal period; 

volume; 

marina width; 

phase angle; 

grid length; and 

frequency 


Subscripts 
high water; 
low water; 
mean tide; 
tidal prism; and 
time step. 


Superscripts 
denotes decreasing values; 
denotes increasing values; and 
denotes constant values. 


MAY 1981 


BREAKING-LIMITED WAVE HEIGHTS 


INTRODUCTION 


Distribution wave heights wave field particular relevance 
ocean engineers. For sea states which the free surface approximately linear 
Gaussian with narrow-band spectrum, this distribution given theoretically 
the Rayleigh law. Numerous field data confirm the validity the Rayleigh 
theory, particularly for low and medium wave heights. However, systematic 
discrepancy between empirical data and the theory often noted towards the 
high wave tail, the theory overpredicting the observations. Likely causes 
this discrepancy include the nonlinear non-Gaussian nature sea waves, 
bandwidth effects, wave breaking, and possibly others. Some studies (3,10) 
demonstrate various compensating procedures for the discrepancy, but fail 
explain why the conventional Rayleigh theory overpredicts field observations. 
Others (6,12,13) examine some effects nonlinearity, finite bandwidth, and 
differences the definitions employed data analysis versus theory. 
particular, Longuet-Higgins (6) and Tayfun (12,13) suggest that the observed 
discrepancy cannot attributed finite-amplitude effect but effects 
associated with finite bandwidth, and differences certain basic definitions 
employed data analysis versus theory. 

appears likely that the skew non-Gaussian nature the sea surface due 
finite-amplitude effects would not directly result reducing theoretical wave 
heights manner consistent with field observations (6,7,12). more plausible 
mechanism wave breaking, which finite-amplitude effect unaccounted 
for the analytical wave models and the development the Rayleigh law. 
previous study the writer (14) qualitatively demonstrated that wave breaking 
will introduce deficiency wave heights towards the high wave extreme. 
This result followed indirect approach which Longuet-Higgins’ (5) approxi- 
mate solution for the joint distribution wave heights and periods was constrained 
satisfy the Stokes limiting steepness deep water. aimed herein 
derive distribution for breaking-limited wave heights following direct 
approach which the statistical description the free surface modified 
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satisfy Stokes-type limiting steepness deep finite water depths. 
the basis this distribution, seen that wave breaking contributes the 
deficiency wave heights towards the high wave tail manner dependent 
various moments the surface spectrum and the still water depth. Comparisons 
with two well-documented field observations are given demonstrate that the 
distribution wave heights limited breaking compares more favorably with 
reality than the conventional Rayleigh law. 


Derinitions AND Resutts 


The oscillations the sea surface given point function time 
represented 


which spectral density the surface; independent and identically 
distributed random phase angles with uniform density function 1/27 [0, 
27]; and radian frequency which related the wave number, 
uniform water depth, through the dispersion relation 


The jih moment the spectrum defined 


0 


spectrum narrow band 


Under this condition, the wave envelope amplitude X(t) defined 
(see, Refs. 12, and 16) 


These definitions suggest that each realization can viewed 
amplitude-modulated wave form with apparent mean frequency i.e.: 


X(t) A(t) cos (wot 


which random phase with the uniform density 1/27 (9). 
Therefore, A(t) and —A(t) represent two symmetrical curves passing through 
wave crests and troughs, respectively; and 


defined the wave height with 


BREAKING-LIMITED WAVE 


recognized the Rayleigh density (5,9,11) with the invariant form 


The preceding formalism underlined with the assumption that the sea surface 
linear and can represented the sum many constituent wavelets with 
different amplitudes, frequencies, and random phases, i.e.: 


N 


The amplitudes are chosen that 


the conditions necessary for the central limit theorem are satisfied, i.e., 
max Aw, with being fixed, and asymptotically approach 
two independent Gaussian variates with zero means and variance, Conse- 
quently, the density both and becomes the Rayleigh law (9,11). Strictly 
speaking, this limiting operation requires that the spectral representation 
and stationary random functions expressed Fourier-Stieltjes integrals 
and respectively. shown Yang (17), the convergence 
and Gaussian variates fairly rapid for small and midrange values 
these variables, but considerably slow towards the extremes. general, 
then, and are not exactly Gaussian unless and the distribution 
wave heights for fixed must deviate from the Rayleigh law particularly 
towards the extreme tail. 

The question which arises now whether there any need all consider 
and opposed their respective theoretical limits, i.e., Eqs. 
and which imply This answered readily noting that 
and require that the range and unrestricted, which not physically 
acceptable. reality, extremes both and are limited whitecapping 
breaking. Therefore, Eqs. and with bounded range should offer 
realistic stochastic representation for the sea surface when constrained 
properly physical grounds. 


deep-water progressive wave, whose limiting form the crest the 
Stokes 120° angle, the acceleration equal (1/2) Thus, typical sea 


Finally, the probability density the scaled heights 
(10) 
N 
n=l 
n=l 
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state whitecaps appear whenever the vertical acceleration the crest approaches 
—(1/2) The limiting form finite water depth given Miche’s (8) 
approximate analysis 


Evidently, Eq. provides semitheoretical transition from the Stokes limiting 
steepness deep water, i.e., 27/7, transitional water depths 
where its validity has been established reasonably well-based field and 
laboratory data (see, e.g., Ref. 15). 

Assume for the moment that the wave described the linear theory. 
this basis, the limiting steepness deep water corresponding downward 


(KH) l 


compared the theoretical value 27/7. Therefore although might 
thought unacceptable use the linear theory for obviously nonlinear effect, 
doing involves error only 11% overestimate. 

The extension these arguments irregular sea state not immediate 
general. However, for sea state characterized with narrow-band spectrum 
such that strictly, each realization the surface has the form 
progressive wave train with slowly varying amplitude, and apparent 
mean frequency, w,. Under these conditions, Eq. should reasonably well 
describe the limiting form irregular wave packet locally. Therefore, 
with being correct second order (12), Eq. becomes 


(KA ) max ko A max yr 


which represents the mean apparent wave number, which related 
through Eq. 


Nonlinear corrections are and, therefore, can neglected 
(12). Now, without any loss generality, let 


implied that the more peaked the spectrum S(w), w,, the smaller 
should used and vice versa. this basis, Eqs. 12, and can 
rewritten 


N 
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which €,; and the scaled wave height when fixed. 
can combined show that 


Thus, inversely related and also depends the relative 
can regarded measure sea state with typical range varying from 
0.1-0.2 for larger waves hurricane much smaller values representative 
gentle seas with apparent whitecaps breakers. Since function 
both and k,D, the distribution narrow-band wave field 
is, essence, dependent the first two spectral moments and the local water 
depth. 


N'? 


Consider now the probability structure can shown (see, e.g., Ref. 
17) that the characteristic function and given 


which and J,(-) the zero-order Bessel function. 
joint density and the inverse Fourier transform Eq. 25, 


which |x| and and the simplification leading the last step 
follows from change variables from (u, polar coordinates. The cumulative 


Using Eq. and paying careful attention the region integration implied 
Eq. 27, can verified that 
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x 


The probability density follows from Eq. differentiation with respect 


This can simplified further, using Heaviside unit step function 
and change variables the integrals, 


VN 
x 


u =—(y2 4 2 


which can substituted Eq. show that 


The exceedance probability distribution 


ww2 
a/2 
cos 
sin 
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The numerical integration Eq. straight forward using Simpson’s 
rule. Fig. shows some typical results constructed this manner for 
10, 15, 30, and 60. The Rayleigh law also included the same figure 
for comparison. observed that the Rayleigh law overpredicts waves with 
heights towards the extreme tail manner increasing with smaller values. 


--- RAYLEIGH 
BREAKING-LIM. 


Wy 


FIG. 1.—Theoretical Exceedance Distribution Breaking-Limited Wave Heights under 
Narrow-Band Conditions and for Different Sea States Implied 


The convergence the Rayleigh limit becomes increasingly slow for larger 
values. 


Comparison Data 


Deep Water Case.—The data examined for this case consist approx 400 
waves collected the Gulf Mexico under reasonably homogeneous conditions. 
contained records Ref. with all pertinent information 
listed Tables and Fig. 5.1 the same reference. These data are considered 


-2 
10 \- 
N=15 
\ N=30 
10° 
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also Longuet-Higgins (5). The empirical exceedance distribution shown 
Fig. was constructed the basis the scatter diagram given Fig. 5.1 
Ref. manner consistent with the notation here. Statistics listed 
(5) 0.234. Therefore, 0.142, and Eq. with k,D requires 
that The corresponding theoretical exceedance distribution given 
Fig. together with the Rayleigh counterpart. evident that the empirical 


DATA 
RAYLEIGH 


BREAKING-LIM. 


FIG. Exceedance Distribution Based Bretschneider Data (1) 
Comparison with Present Theory and Rayleigh Distribution 


points are approximated the distribution more favorably than the 
conventional Rayleigh law. Furthermore, note that 3.162 
theoretically, comparison with the observed 2.198. 

Finite Water Depth Case.—The data for this case are from Dattari, et. 
Unfortunately, the data represent collective summary different records, 
each containing approx 100 waves. Statistics pertinent each individual record 
are summarized Table Ref. which indicates variability between 


) 

‘ 

10° 

\ 
\ 
\ 
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for the majority the records. Therefore, 0.036, and Eq. 
requires that 30. The empirical exceedance distribution constructed 
the basis Fig. Ref. shown Fig. here together with the distributions 
corresponding and the Rayleigh law. The comparison indicates again 
that the distribution with more favorable than the conventional Rayleigh 
This cannot considered too contradictory considering that the actual sample 
size approx 100 waves for each constituent record. may interest 


DATA 
--- RAYLEIGH 


FIG. Exceedance Distribution Based Data Dattari et. (2) 
Comparison with Present Theory 30) and Rayleigh Distribution 


also indicate that these waves were deep water Eq. would require 
78, and the corresponding distribution would shifted slightly towards 
the Rayleigh limit. 

The wave heights these two cases and almost all other field observations 
are based conventionally zero-crossing definition. The Rayleigh law and 
the distribution derived herein are relevant envelope amplitude definition. 
The two have fundamental difference. Specifically, zero-crossing wave height 
represent two points &(t) separated half the mean period, 


q 
\ 
\ 
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such this definition represents averaging operation damping variations 
about its mean. demonstrated the writer (13), the resulting sample 
space must have the same mean, show excess waves with heights about 
the mean, and deficiency towards both extremes away from it. Therefore, 
histograms and empirical exceedance distributions constructed from zero-crossing 
wave heights could display the discrepancy question manner further 
enhanced this difference definitions. fact, this appears quite 
likely for the aforementioned deep water case because the corresponding 
wave height histogram given Fig. Ref. displays these characteristics. 


Summary 


Empirical distributions derived from field data differ from the theoretical 
distribution wave heights given the Rayleigh law. principal cause 
this discrepancy suggested that the conventional theory does not 
any way account for wave breaking. Based the assumption that the sea 
surface characterized with narrow-band spectrum, was shown that the 
distribution wave heights limited breaking provides probability description 
which compares with two well-documented field observations more favorably 
than the Rayleigh law. 

The breaking-limited distribution depends the first two spectral moments 
the sea surface and the local water depth, and essence, represents 
generalization the conventional theory. Therefore, for certain combination 
these parameters, for larger values representative rough 
sea states and shallower depths where wave breaking likely observed, 
can deviate from the Rayleigh theory substantially. the other hand, for 
small values typical gentle seas with apparent breakers, 
becomes identical with this conventional theory. 

Whether these comparisons and conclusions can generalized with confidence 
remains seen through comparisons the proposed theory with additional 
field lab data, both, reliable quality. ths respect, imperative 
ensure that empirical data are sufficient size and consistent with the definition 
wave heights based the wave envelope. 
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INTRODUCTION 


This paper concerns 4-year experiment Hawaiian waters dealing with 
the forces exerted long-period swell test pipe mounted the sea 
floor (11.3 water. Four separate angles were involved between 
the wave fronts and the pipe: 0°, 15°, 35°, and 52°. The wave-induced water 
motion near the test structure was measured were the horizontal and vertical 
forces instrumented portion the test pipe. Although two earlier 
presentations (7,8) have dealt with data obtained this investigation, has 
not been possible consider the complete set data until this point. Instanta- 
neous force values times zero velocity have been used with concurrent 
accelerations derive information (horizontal) inertia coefficients. These 
correspond closely ideal flow predictions. The focus this paper 
peak forces, however, and maximum force coefficients have been obtained 
normalizing peak forces for drag coefficients. The trends these horizontal 
and vertical maximum force coefficients with changes wave-pipe angle and 
(adapted) period parameter are examined for the data both mean-value and 
extreme-value forms. Results from the experiment are used simulated 
submarine pipeline design problem. 


AND 


Introduction.—Consider rigid pipe outside diameter set parallel and 
distance above level sea floor. the context this paper, the intent 
determine the peak horizontal and vertical forces short length 
this pipe during the passage certain surface wave defined its height 
its period the water depth the location the pipe, and the angle 
the wave front with the pipe, The maximum horizontal force, perpendicular 
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the pipe, that occurs during the passage this wave will represented 
the peak vertical force will written 

The maximum water velocity and acceleration near the sea floor could 
predicted for the wave described above through use any least dozen 
wave theories (11). Let the theoretical predictions the peak velocity and 
acceleration any theory represented and respectively; the 
symbols for the real maximum velocity and acceleration will and 
Kinematics used this paper will only the latter type, making the assumption 
that measured kinematics are fact the actual ones. Every wave theory errs 
least some degree predicting and (1), and correction factors 
are necessary for and 

Table contains relevant information obtained the end this project 
(12). The data were obtained for the following wave conditions: 3.0 
8.2 (0.9-2.5 m); 5.0 (1.5 m); 11.5 17.4 sec; 16.0 sec. 
(Overbars refer sample mean values.) The kinematical data can described 
by: 1.19 3.57 fps (0.36-1.09 m/sec); 2.33 fps (0.71 m/sec); 


0.27 1.97 ft/sec? (0.08-0.60 0.95 ft/sec? (0.29 
The symbol Table represents the standard error estimate. 


TABLE 1.—Comparison Predicted and Measured (Ducted Meter) Peak Kinematics 


Prediction Prediction the 
Kinematical Airy Theory Stream Function Theory 


quantity 
(1) 
6.6% low 0.31 fps 10.7% high 0.31 fps 
(0.095 m/sec) (0.095 m/sec) 
11.2% low 0.26 4.3% high 0.25 
(0.079 (0.076 


Peak Forces.—The equations that will used connect the peak kinematics, 
pipe characteristics, and the maximum forces are the following (4): 


which the water density. Formulations other than Eqs. and are 
possible (2), but the writers have found these the most useful for the 
data herein. apparent that these two equations contain explicit mention 
the peak acceleration but there implicit inclusion this variable 
through the realization that the maximum force coefficients and 
must functions the relative effect acceleration. The parameter that 


will used indicate such tendency will adapted period parameter 
given the equation 
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Parameters and Instantaneous Inertia Force parameters 
interest are follows: 


max 


the Reynolds number 
v 
which the kinematic viscosity the water; 


the relative pipe clearance 


the relative roughness the pipe surface 


which the average height surface protuberances; 


the relative distance water particle motion 


which ideally the horizontal distance water particle movement. 
estimate was obtained this investigation graphically integrating under 


3.4 


maximum 3.29 
minimum 2.00 


0.2 
A/D 


FIG. Flow Inertia Coefficient for Pipe Function Pipe Clearance 
from the Boundary and Angle Attack 


€ 
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U-time curve from the instant when the time interest, e.g., 


max* 


The instantaneous inertia coefficient 


which the flow acceleration when The ideal-flow counterpart 


will written and this coefficient function both A/D and 


Experimental in-depth account the design, fabrication, 
installation, operation, maintenance, and diving aspects this project presented 


TABLE Test Conditions 


degrees 
(4) 

3,5,9 

16,19,20 

22,23 

29,30,31,32,33,34 

39,40,42,43,44 

45,46 


Grace (5,6). The following summary that material. The test site 
was (11.3 water offshore from Honolulu and the sea floor 
the area the test site was fairly flat coral rock. The tidal range Hawaii 
very small, averaging 1.2 (0.4 m); near-bottom tidal currents the site 
were negligible. The trade wind fetch the site was only 2,000 (650 m). 

test pipe structure was designed and fabricated, then installed thin 
gravel bed placed the natural sea floor even out the subtle peaks and 
hollows. The test pipe structure consisted two parts—first, 16-in. diam 
(406-mm) steel pipe, 17.5 (5.3 long, and second, steel and concrete 
base which the pipe was attached. The base measured in. 
(5.0 2.4 plan and was faired with the natural sea floor through 
use dozens concrete beams, quick-setting underwater cement and epoxy, 
and steel chain. Before the test pipe system was placed the ocean, the direction 
propagation Southern Hemisphere swell the site was established (+5°) 
through extensive sea floor observation water motion. 

Only 40-in.-long (1.04-m) near-central portion the pipe was used 


No. 
waves 
Surface pro- 
Phase A/D finish cessed 
(1) (5) (6) (7) (8) 
0.031 ribs 0.016 219 
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force sensor. This segment was attached the two side portions the pipe, 
welded the base, through two internal strain gage beams. One these reacted 
vertical forces and the other horizontal ones. test day the vertical 
channel was calibrated placing known weights the test section and the 
horizontal channel was calibrated using shunt resistor coupled with the 
vertical force calibration. The latter technique had been checked earlier the 
field, exerting known horizontal force the test section, and found 
accurate within 2%. 

Several different pipe arrangements were used during the years experi- 
mentation. the first year, the pipe was oriented perpendicular the wave-in- 
duced water motion; the pipe clearance from its base was in. (76 mm). The 
pipe base was 4.5 in. (114 mm) deep. Throught the remainder the project, 


pulse history from ducted meter [nominal calibration is 0.563 ft 
(0.172 m) of water motion between pulses) 


force \offshore 


horizontal 
force 


sec 


crest 39-850 
trough 39-85! 


crest 39-852 
trough 39-853 


force up 


FIG. 2.—Example Strip Chart Record Concurrent Variation Ducted Meter 
Flow Speed Output and Horizontal and Vertical Forces Pipe Test Section 


the pipe clearance was changed 0.5 in. (13 mm). This was done laying 
keyed, concrete slabs over the old base and the base depth thus became 
in. (178 mm). The original 3-in. (76-mm) clearance, found unnecessary, 
was used ensure that coral rock rubble, conceivably swept over the sea 
floor during heavy wave action, would not jam the test section. The 0.5-in. 
(13-mm) gap was used the best compromise between having the pipe directly 
the base, truly simulate real pipe, and have sufficiently far from 
the base for the test section function without any possibility 
Observations sand near the test section provided strong evidence that there 
was perceptible water motion through the 0.5-in. (13-mm) opening. 

Initially, the pipe was used condition. Midway through the 
second year work, artificial roughness elements (Table were added 
the pipe test section and adjacent parts the flanking pipes, and these 
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were maintained through the rest the project. The angle between the nominal 
wave fronts and the pipe during the second year remained 0°, but for the 
third year angle 15° was employed. Two separate angles were used during 
the final year work, first 35° and then 52°. The pipe reorientations were 


TABLE 3.—Peak Kinematical and Force Values for Wave Train Fig. 


feet 
per 
second 
(meters 
Wave per 
feature second) 


(1) 


degrees 


(2) 


feet 
per 
second 
squared 
(meters 
per 
second 
squared) 


Sample 
size 


(3) 


feet 
per 
second 
(meters 
per pounds 
second) (newtons) 


(5) 


Standard 
deviation 


(5) 


pounds 
(newtons) 


(6) 


Maximum 


(7) 


carried out, after removal the fairing beams, through use large, air-filled 
lift bag and pulling from one the four anchor blocks 
the sea floor used keep the project boat anchored over the test site. 
special vane was used measure the angle between the pipe and the direction 


(2) 
39-850 1.49 0.85 1.91 31.6 26.0 
(0.45) (0.26) (0.58) (142) (116) 
39-851 1.91 0.91 1.75 35.6 48.8 
(0.58) (0.28) (0.53) (159) (219) 
39-852 1.75 1.10 1.99 48.8 48.0 
(0.53) (0.34) (0.61) (219) (215) 
39-853 1.99 0.70 1.40 22.4 18.4 
(0.61) (0.21) (0.43) (100) (82) 
TABLE 4.—Peak Force Data Summary for Phase 
Maximum Force Data, pounds (newtons) 
Force Mean 
(1) (4) 
212 60.9 25.5 19.5 148.8 
(271) (113) (87) (662) 
210 57.4 21.4 18.5 128.1 
(255) (95) (82) (570) 
243 38.3 16.1 13.2 116.6 
(170) (72) (59) (519) 
229 32.4 9.6 12.2 59.0 
(144) (43) (54) (263) 
195 24.0 8.1 10.3 57.7 
(107) (36) (46) (257) 
186 18.4 6.5 6.9 40.0 
(82) (29) (31) (178) 
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wave-induced water motion any Data Day. 

During data-taking days (see Table 2), data were recorded the project 
vessel its forward cabin. Electrical cables connected strip chart recorders, 
powered portable alternator, with the sensors the sea floor. Two channels 
data involved the horizontal and vertical forces the pipe test section, 
coming from the strain gage half-bridges, whereas third channel involved 
the output from the ducted current meter. This record involved series 
pulses, magnets the tips the three blades the ducted meter impeller 
spun past enclosed reed switch the side the meter. All three time 
histories were recorded the same chart, shown Fig. retraced raw 
data history for Data Day 39. Table provides some details the waves 


TABLE 5.—Further Information for Phase Data 


degrees 


pounds 
(newtons) (newtons) 


(4) (5) 


shown Fig. this table represents the peak flow speed under the 
preceding wave feature (i.e., crest trough). 

Force Data.—The ranges maximum horizontal and vertical forces measured 
during Phase (Table are presented Table Additional information 
general nature shown Table that table, the following symbols are 
used: sample size; correlation coefficient between same-wave pairs 
the time delay relative the same-wave and s.d. represents 


max? 


standard deviation. Note that negative values Table indicate leading 


max 


The time delay with respect was also investigated for the 
Data Day data from Phase 0°). sample size produced 


seconds 
Mean s.d. mum 
(1) (6) (7) (9) 
203 0.90 0.89 210 0.45 0.27 —0.40 3.05 
174 0.61 0.90 164 0.49 0.37 —1.28 1.71 
181 0.65 0.74 176 0.48 0.49 3.17 
TABLE 6.—Correlation Peak Horizontal and Vertical Pipe Forces for Phase 
Test 
sequence 
(1) 
106 0.50 35.8 21.0 
(159) (93) 
0.80 58.2 55.2 
(259) (246) 
0.94 68.5 66.0 
(305) (294) 
P’ 
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mean 0.45 sec, standard deviation 0.45 sec, and range from —0.60 
2.90 sec. There were only three negative lags and the high lag 2.90 sec 
was one only two time intervals, both for troughs, larger than 0.90 sec. 
Fortyone (75%) the Data Day time delays fell the interval 0.20 
0.50 sec. 

The correlation between the peak horizontal and vertical force data from 
Phase presented Table The results suggest that increasing pipe roughness 
and decreasing pipe clearance promote coupling between maximum horizontal 
and vertical wave-induced forces for pipe near boundary. 


Errects 


major aim this research was chart the manner which the force 
coefficients (Eq. and (Eq. vary with A/D, and (factors 


TABLE 7.—Pipe Force Differences Due Factors Other Than Gross Kinematics Under 
Particular Wave 


feet 
per 
feet second feet 
per squared per 


second per pounds pounds 
no. second) squared) second) seconds (newtons) (newtons) 

(2) (3) (5) (6) (7) (8) 


under control) well Reynolds number and adapted period 
(parameters determined the wave environment). stressed, however, 
that even for the same test day (same A/D, and, within limits, and 
for two waves having the same and (thus, same Reynolds number 
and period parameter say nothing the same squared velocity appearing 
and 2), forces are different. 

The fact that flow history appears play strong role such variability 
from zero peak flow speed. Although some the scatter the force coefficient 
data, shown presently, may due experimental error, the writers 
attribute the bulk such variability effects beyond the factors included 
mathematically the analysis. particular, flow history and flow turbulence, 


520 2.55 1.36 2.98 2.21 76.1 
(0.78) (0.41) (0.91) (341) (357) 
770 1.47 1.31 2.93 2.14 60.7 58.0 
(0.45) (0.40) (0.89) (272) (260) 
044 2.67 1.41 2.24 1.65 74.5 
(0.81) (0.43) (0.68) (334) (334) 
050 2.10 1.39 2.24 2.47 62.0 62.3 
(0.64) (0.42) (0.68) (278) (279) 
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probably intimately related, play strong dispersive roles. Because such random 
effects, force coefficients will considered later two senses—mean trends 
and extreme values. 


Table summarizes the information obtained the project for the instanta- 
neous inertia coefficient given Eq. The distributions values for 
the data sequences shown Table appeared adequately fitted the normal 


TABLE and Instantaneous Inertia Coefficients 
Instantaneous 
Inertia Coefficients Correlation 


Test 95% confidence Coefficients 
sequence Mean interval With 
(4) (6) (7) 


1.98-2.26 —0.23 
2.27-2.55 0.28 
2.50-2.90 0.08 
2.63-2.85 —0.09 
2.31-2.45 —0.04 
1.85-2.05 


distribution. The mean coefficient values come satisfactorily close the 
theoretical values, (Fig. 1). 


Forces 


Introduction.— Values plotted versus adapted period parameter, 
for example Figs. 3-6, yield clouds points. Two approaches will 
used later with these data: the one hand, mean line will derived; 
the other, envelope trend will considered. Although some work was done 
with the distribution residuals about the mean-line curves, this approach 
will not used herein the dispersion the distribution was correlated 
with the independent variable (in predictor sense) 

The correlation coefficient between and Test Sequence 2A, for example, 
(211 points) —0.75. But the correlation coefficient between and 
for these same data —0.68. then logical inquire whether 
really varying negative way with both these parameters perhaps 
varying such way with only one and the two parameters are themselves 
positively correlated. For these data the correlation coefficient between and 
0.85. 

answer this question, (also was plotted both against for 
specific ranges and again against for specific ranges. Adapted period 
parameter classes were taken 0.2-unit width with ranges 1.0-2.5, 
2.5-4.0, and >4.0. The data showed clearly that there were visible trends 
purely the basis thus the data have been presented Figs. 3-6, 


2.65 218 0.81 0.06 
1.69 195 0.71 0.18 
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i.e., versus identification the individual points the applicable 
range N,. similar approach, but statistical, appears Tables and 10. 
The case, missing these tables, discussed Appendix Table 
presents general correlation information for Phase data. 


FIG. 3.—Variation Maximum Horizontal Pipe Force Coefficient with Adapted Period 
Parameter 35°) 


FIG. 4.—Variation Maximum Horizontal Pipe Force Coefficient with Adapted Period 
52°) 


Maximum Force Coefficients: Mean Trends.—The following forms expres- 
sions have been found useful fitting the classified mean trends for 
and versus adapted period parameter: 
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e ar 56:2 


FIG. 5.—Variation Maximum Vertical Pipe Force Coefficient with Adapted Period 
Parameter 15°) 


FIG. 6.—Variation Maximum Vertical Pipe Force Coefficient with Adapted Period 
Parameter 35°) 


Gam = Cyt 


The velocity-dependent horizontal force coefficient Eqs. and came 
from analysis wherein peak forces were predicted Morison form 
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TABLE 9.—Mean and Extreme Values 


Mean Extreme Mean Extreme 
(3) (4) (7) (9) (10) 


4 
4 


N 


ww2 
3.2 2.48 3.39 2.05 2.81 1.43 
5.0 2.06 2.29 1.69 1.03 
5.6 1.60 1.88 1.39 0.86 
5.8 2.01 2.01 0.65 0.75 
6.0 1.84 2.01 1.28 0.90 1.32 
6.6 1.63 1.88 1.61 0.50 0.50 
8.8 
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TABLE 10.—Mean and Extreme Values 


Mean 
(6) (7) 


ww2 
(1) (4) (9) (10) 
0.8 9.32 9.32 7.53 8.93 3.18 
1.2 3.69 6.09 5.80 2.33 4.02 
1.6 3.80 4.89 4.06 5.44 1.88 4.39 
1.8 3.69 4.63 2.87 5.14 1.60 2.28 
2.0 3.58 5.41 4.15 1.32 3.15 
2.2 3.50 4.65 2.43 3.33 1.71 2.52 
2.6 2.81 3.63 2.13 3.22 3.07 
2.8 2.95 4.09 2.02 3.23 1.21 1.98 
3.4 2.00 2.62 1.71 2.39 0.94 1.58 
3.6 2.93 1.35 1.86 1.60 
3.8 2.16 2.67 1.56 0.80 1.16 
4.4 1.79 2.41 1.19 1.77 0.98 1.12 
4.6 1.97 2.57 1.28 1.60 0.73 1.13 
4.8 1.94 2.64 1.26 1.66 0.75 1.01 
5.4 1.06 1.19 0.87 1.49 
5.6 1.58 0.92 1.21 0.65 0.85 
5.8 1.87 1.87 1.06 1.39 0.57 0.66 
6.0 1.48 1.59 1.10 1.10 0.76 1.31 
6.4 1.25 1.46 0.84 1.04 
6.6 1.57 0.98 1.27 0.50 0.50 
7.0 1.29 1.29 0.81 0.81 
7.2 1.29 1.29 0.99 0.99 
8.8 218 200 
7 
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the equation where was the acceleration-related coefficient and was 
the velocity-dependent one. The use this coefficient the asymptote for 
the maximum horizontal force coefficient certainly reasonable. was used 
asymptote for the maximum vertical force coefficient since appeared, 
after various trials using other values, the best choice. 

The coefficients, and appearing Eqs. and were chosen from 
the three horizontal and vertical data sets Phase they are also listed 
Table 12. The fit the resulting six lines the mean trends the 
and data was excellent, and Figs. and present the resulting curves. 
Two comments are perhaps order assist understanding these two figures. 
the one hand, the same data classes were used for Tables and 10, 


TABLE 11.—Correlation Information for Phase Data 


Coefficient 
Correlation Coefficient with Parameter 


degrees parameter 
(1) 


TABLE for Eqs. and 


Test 
sequence 


0.6, 0.6 0.8, etc. Secondly, the open symbols the 
figures indicate that only one value was involved these classes; thus, the 
particular values not have the same weight those obtained other classes 
from samples size two greater. 

interest compare the and values for 35° and 
52° those obtained for 15° or, turns out, those 
values applicable for 0°. all four cases, there apparent trend 
the ratios with increasing adapted period parameter. Summary results for 
the four cases are presented Table 13. Note that situations involving only 
one value class were ignored. 

Maximum Force Coefficients: Extreme Values.—It also interest consider 


(4) (5) 
210 
219 0.85 1.00 
298 0.84 1.00 
200 0.84 1.00 
Ch A ! G, A 2 G, 
(1) (2) (4) (5) (7) (8) 
2.87 0.239 1.99 6.50 0.478 0.38 
0.77 2.20 0.222 1.81 4.95 0.525 1.39 
0.41 2.26 0.285 1.16 2.03 0.337 0.82 
) 
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the extreme values the maximum force coefficients Tables and 10. 
For each class each table, excepting those situations where there only 
one point, the ratios obtained dividing the extreme coefficient value 
the mean have been computed. Down playing the tendency towards lower ratios 
for high values, because the small class sample sizes, there appear 


Dete off Chart 
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FIG. 7.—Variation Class-Average Maximum Horizontal Pipe Force Coefficient with 
Adapted Period Parameter for Three Test Pipe Orientations 


Doto off Chort ¥-08-10 


FIG. 8.—Variation Class-Average Maximum Vertical Pipe Force Coefficient with 
Adapted Period Parameter for Three Test Pipe Orientations 


tendencies among the six separate groups towards systematic down 
trends with increasing values. For that reason, convement view 
overall information and the results are given Table 14. clear that for 
both horizontal and vertical maximum force coefficients, the averages increase 
the angle increases. 
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The Question Very complete for use 
computing peak horizontal forces are the following: 


should stressed that work done this project concerns the latter case; 
the former situation theoretical origin any event. For short region 


TABLE 13.—Ratios Mean Maximum Force Coefficients Within Classes 


degrees Mean 
(4) 


Ratio 


max extreme 


15°) 0.63 


max 


The point intersection between the hyperbolic relationship given Eq. 
with the mean lines shown Fig. theoretically gives estimate the cutoff 
value, between the two forms Eq. 11. Unfortunately, for each Test 
Sequences 2A, 2B, and 2C, the two lines not actually cross for low values 
but virtually parallel each other. terms both the closeness the lines 
and the paucity points this investigation for proposed that 
1.0. Thus, Eq. can rewritten 


Ratio Low High 
(1) (5) (6) 
TABLE 14.—Ratios Extreme Mean Maximum Force Coefficients 
Low High 
(6) (7) 
0.10 1.42 
1.30 0.14 1.06 1.65 
1.40 0.16 1.04 1.68 
1.29 0.13 1.07 1.65 
1.44 0.19 1.12 1.79 
1.61 0.33 1.14 2.39 
point overlap between the two relationships, can shown that 
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Note that there relationship such the foregoing for the vertical maximum 
force coefficient. 


Introduction.—One option considered for the sewer outfall Wanganui, New 
Zealand, was post-tensioned concrete pipe, in. (1.7 outside diameter. 
possibility this case was that the line laid along the seabed some 6200 
(1900 out water depth (12 m). The pipe would stapled 
the seabed required. Storm conditions the Tasman Sea, which offlies 
Wanganui, are such that was reasonable assess tha planned outfall’s stability 
breaking waves. The steepest angle approaching wave orthogonals the 
pipe, 35°, was obtained from refraction analysis involving waves approaching 
from due west. Thus, for design, 55°. 

the present section, the results this research investigation Hawaii 
will used evaluate the peak horizontal and vertical forces the (assumed 
rigid) Wanganui outfall (12 water. Some consideration will also 
given the distribution force along the outfall this depth. The problem 
computing forces under waves that are the process breaking rather 
than incipiently breaking has been described Grace (4) and will not covered 
here. 

Kinematics, 10-sec period for the waves incipiently 
breaking approximately (12 water towards the end the outfall. 
Although has been shown (1,3,9,11,13) that Airy theory accurately predicts 
near-bottom velocities under the crests even near-breaking waves, the 
prediction the trough velocity too high and the peak acceleration too low. 
Use tabulated (1) stream function theory information and extensive interpola- 
tion gives 
—4.84 fps (—1.48 m/sec) 
5.48 ft/sec? (1.67 m/sec’) 


max 8.67 fps (2.64 m/sec) 


which and the trough and crest velocities, respectively, and 
the peak water particle acceleration occurring between trough and crest 
between crest and trough. Assume, for simplicity, that these three theoretical 
kinematical values are equal the real quantities and 


max? 


respectively. 

The values adapted period parameter for the trough and crest are, respec- 
tively, 


2 
max 


and 


max 


Both valves Eq. are necessary since the distribution force the 
pipe should computed. The value indicates use Eq. 13a, whereas 
the value points Eq. For the trough case, Fig. (A/D 0.00 
and 55°) yields 1.9, and Eq. 13a gives 


which the presubscript indicates trough. the crest case, interpolation 


and extrapolation Figs. and yield 1.7 and 1.3, respectively 
(for 2.46). Thus, 


max 


, Pp 2 


Since formulation such Eq. 13a available for vertical forces, the 
assumption will made that 


, 
t max 


For pipes very close proximity the bottom, the approximate validity 
Eq. 19a was assured information gathered during this project. 

The wave length (L) the 30.4-ft (9.3-m) high design wave is, according 
the stream function theory, 344 (105 m). the depth were reasonably 
constant, the wave feature would theoretically repeat itself intervals 


along the pipe. The depth change over this distance approximately 
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that the original assumption constant depth adequately substantiated. 
Assume wave crest the offshore end the pipe. Horizontal forces 
710 (10.3 and vertical forces 540 (7.9 kN/m) exist 
the end the pipe and 420 (128 inshore that location; horizontal 
and vertical forces /m) and (7.7kN/m), respectively, 
occur 210 (64 shoreward the end the outfall. The horizontal crest 
and trough forces act opposite directions and the horizontal force distribution 
along the pipe can approximated fairing smooth curve that joins 
these points. The vertical crest and trough forces both act upward. The distribution 
vertical force along the pipe passes through these values and dips between. 
This minimum vertical force can assumed equal zero shown Fig. 


Ocean experiments with test pipe have resulted extensive wave force 
data for angles between wave fronts and the pipe 0°, 15°, 35°, and 52°. 
The horizontal and vertical pipe force data have been converted into maximum 
force coefficients through use measured kinematics and these coefficients 
have been related values adapted period parameter. This information 
useful engineers appraising peak wave forces exposed and rigid submarine 
pipes. 
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This appendix included order explain the absence the Phase 
0°) force coefficient data from Tables and and Figs. and 
has been shown Table that the three parts Phase had different 
combinations. However, the force coefficient data from 
Test Sequences and (different roughness only) were indistinguishable 
and the two sets data have been pooled what follows. The combined 
and information [pipe clearance 0.5 in. (13 mm)] reproduced 
the format Tables and Table 15. Comparison these data with 
those for 15° Tables and shows clearly that the trends Table 
lie somewhat lower than those Test Sequence 2A. the average (excepting 
classes with the coefficient ratios are 0.83 and 0.85 the horizontal 
and vertical directions, respectively. 

The finding that the ratio these coefficients not very close to, and 


| 
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TABLE 15.—Classified Mean and Extreme Behavior Maximum Force Coefficient 
for Test Sequences and 


Extreme 


(7) 


Ae 
Akar 


4 
4 


w 
w 
N 
NR 


(1) (3) (4) (6) 
0.4 
0.6 
0.8 3.34 3.37 
1.0 
1.2 3.67 3.67 
3.03 3.86 
1.6 3.02 3.68 
1.8 3.21 4.14 
2.0 2.70 3.58 2.70 4.00 
2.73 3.74 3.00 4.33 
2.4 2.37 2.68 3.04 3.26 
2.6 2.35 3.45 2.43 4.00 
2.8 2.16 2.73 2.25 3.33 
3.0 2.63 3.14 2.37 2.97 
3.2 2.11 2.86 1.86 3.02 
3.4 2.24 2.99 2.38 3.21 
3.6 1.78 2.33 1.54 2.43 
3.8 1.95 2.39 1.55 1.99 
4.0 1.89 2.53 1.29 1.65 
2.04 2.52 1.78 1.87 
1.81 2.21 1.63 2.2 
2.10 2.89 1.69 2.14 
1.73 1.97 1.42 2.04 
5.2 1.59 2.21 1.30 1.91 
5.4 1.34 1.50 1.07 1.13 
5.6 1.72 2.00 1.87 2.32 
5.8 1.69 1.70 1.57 1.70 
6.0 1.19 1.35 1.04 1.21 
6.4 1.40 1.65 1.10 1.22 
6.6 1.42 1.42 0.79 0.79 
6.8 1.64 1.64 1.46 1.46 
7.0 1.29 1.36 0.99 1.07 
8.2 1.16 1.16 0.80 0.80 
8.8 
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perhaps somewhat greater than, unity surprising. However, the impression 
the writers that the and 15° cases are essentially the same and 
that the 0.83 and 0.85 ratios resulted from irregularity the data gathering 
and reduction. The ducted current meter used this project has different 
calibration (10%) depending upon whether flow enters one end the other 
the duct. That this difference was large 10% did not surface until 
detailed unsteady-flow calibration the meter (which, incidentally, demonstrated 
the excellent dynamic response the sensor) was carried out between Phases 
and Whereas very careful attention the fore and aft orientation the 
meter was maintained throughout Phase this was unfortunately not always 
the case during Phase The writers believe that during one and perhaps two 
the Data Days Test Sequences and Table the meter orientation 
was inadvertently switched. This would have resulted flow speeds too high 
for the crest case, the primary one for which data were processed, leading 
coefficients that were too low. 

Before leaving the Phase data, worthwhile considering differences among 
the 1A, 1B, and data sets. Table presents the different conditions for 
the three parts Phase 

There were discernable mean-trend extreme-value differences for the 
horizontal maximum force coefficients across Test Sequences 1A, 1B, and IC. 
There were detectable extreme-value differences across the three data sets 
for vertical peak force coefficients, and the same applied the and 
mean-value trends. However, the mean-value line clearly lay below those 
for Test sequences and IC. 

These overall results indicate, first, the very small role that reasonable-scale 
roughness plays determining wave forces real ocean submarine pipes. 
The writers believe that the heavy wake and nonwake turbulence swept back 
and forth over uneven sea floor responsible for this result. other words, 
turbulence effects overwhelm roughness effects. The second observation con- 
cerns pipe clearance. Although the lack change horizontal peak force 
coefficients might have been anticipated from steady flow results (10), the absence 
change the extreme values for the vertical maximum force coefficients 
unexpected and important result. 
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LONG-TERM DISTRIBUTIONS OCEAN WAVES: 
REVIEW 


Michael St. ASCE and Neil MacKenzie” 


INTRODUCTION 


Important advances have now been made into various aspects the design 
offshore structures with respect wave action. These include: (1) The 
development highly accurate wave theories; (2) the prediction wave forces 
slender structural elements terms the incident fluid flow, and large 
bodies the basis diffraction theory; (3) realistic wave hindcasting from 
known wind field data; and (4) spectral descriptions wave and wave loading 
variations during given sea state. One important step often encountered 
design the estimation extreme design wave (or sea state) the basis 
recorded hindcast wave data. This generally involves selecting and fitting 
suitable probability distribution wave height data, and extrapolating this 
locate suitable design value. However, even though the procedure 
considerable practical importance, there appears systematic account the 
literature the various steps involved applied all the different distributions 
that have generally been adopted. 

The selection design wave from series wave records usually 
carried out the following stages: 


Data consisting wave heights and periods are collected over long 
time (e.g., years) the site interest. Alternatively, hindcasting technique 
may used provide wave height data over much longer time span (say 
more). 

plotting formula used reduce the data set points describing 
the long-term probability distribution wave heights. 

These points are plotted extreme value probability paper corresponding 
chosen probability distribution function. 

straight line fitted, often eye, through the points represent 
trend. 


Prof., Dept. Civ. Engrg., Univ. British Columbia, Vancouver, British 
Columbia, V6T Canada. 

Engr., Swan Wooster Engrg. Co., Vancouver, British Columbia, Canada. 

Note.—Discussion open until October 1981. extend the closing date one month, 
written request must filed with the Manager Technical and Professional Publications, 
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This paper part the Journal the Waterway, Port, Coastal and Ocean Division, 
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The line then extrapolated locate design value corresponding 
chosen return period, T,, chosen encounter probability, 


indication some the preceding steps provided the wave height 
exceedence diagram sketched Fig. (Relationships between alternative 
ordinate scales based the reduced variety, the exceedence, the return 
period, and the encounter probability, are shown. The scale for 
derived for design lifetime, yr.) Much the information needed 
carry out the aforementioned procedure widely scattered the literature, 
often being described any one case for one specific probability distribution 
and fitting technique only. Furthermore, emphasized Nolte (23), quite 
different models have been used basis for defining the individual data 
points which distribution fitted. The present paper serves assemble 
the results needed applying any the more common alternative distributions 


0.001 
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FIG. 1.—Wave-Height Exceedence Diagram Based Log-Normal Distribution 


used, and review the estimation procedure with its principal variations and 
areas difficulty. The emphasis the paper given the selection and 
fitting probability distribution the data. The complimentary step wave 
hindcasting from known wind field not itself treated here, although the 
extrapolation such hindcast data extreme conditions covered part 
the general procedure described. 


Extreme 


Several probability distributions have been used proposed describe extreme 
wave statistics. These include the log-normal distribution, and the Extremal 
Types and III (both lower-bound and upper-bound) probability distributions. 
Although these all have theoretical base, they are used here essentially 
empirical fit the data. 


0.01 

0.10 
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the following, the cumulative probability, P(H) Prob(H(t) H), 
the probability that the argument not exceeded any randomly chosen 
value, H(t), and the probability density, p(H) dP(H)/dH. describing the 
distributions, convenient adopt the following notation for the parameters 
used define any specific distribution: the shape parameter which determines 
the basic shape particular distribution; the scale parameter which 
controls the degree spread along the abscissa (variate axis); and the 
location parameter which locates the position the density function along 
the abscissa. the special case the Type III distributions described, 
locates one end the density function. 

Table gives the expressions for P(H) defining the different distributions 
considered here, and also includes expressions for their means and variances. 
Sketches each the corresponding probability densities, plotted with linear 
scales along both axes, are shown Fig. for various parameter values. 

common practice plot any set measured data such that the selected 
distribution lies straight line, since this aids visualization the extrapolation 
procedure and permits eye fitting where appropriate. Thus, particular probability 
paper first selected, equivalently corresponding scales are constructed, 
order meet this requirement. The linear ordinate scale, related 
the cumulative probability, and the linear abscissa scale, related 
the variate, according the relationships given Table There now 
intercept which are given terms the parameters and indicated 
the table. the light the preceding comments, each distribution now 
briefly mentioned turn. 


TABLE 1.—Probability Distributions Used Describe Long-Term Wave Heights 


Distribution Range T J Cumulative probability. P(H) 1 Mean Variance 
(1) (2) (3) (4) (5) 


Lognormal (1/V 2m) -! exp (20 + a’)fexp @’)- 
o ah 2 a 2 


exo { -( «+70 —e* 
(=« + 0.586) 


(= 1.6407) 


Type Ill, 
(Lower 


Bound) 


Type Ill, 
(Upper 
Bound) 
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Log-Normal Distribution.—The log-normal distribution corresponds 
possessing normal distribution. was the first distribution fitted 
long-term wave height data (see Ref. 19) and has commonly been used extreme 
wave prediction following the work Draper (8). Some examples illustrating 
this approach include those given Draper (9), Powers, Draper, and Briggs 
(28), Singh and Draper (32), Draper and Driver (10), and Dattari (7). These 
examples and others have been compared Ouellet (25). 


FIG. Densities Selected Distributions Plotted with Linear Scales: 
(a) Log-Normal; Type (c) Type Type (e) Type 


Typical densities with linear scaling and with different values are shown 
Fig. 2(a). The corresponding log-normal probability paper readily available, 
and easy construct since the ordinate scale, relating given 
Table listed probability tables (e.g., see Ref. 1). 

Extremal Type Distribution (Also Termed Gumbel Fisher-Tippett Distribu- 
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(15) developed this distribution considerable degree the 
context flood prediction and also frequently used extreme wind prediction 
(e.g., see Ref. 31). The Type distribution has been examined the context 
extreme wave prediction St. Denis (29,30), and example its use 
has been given Khanna and Andru (20). 

sketch the density function with linear scaling shown Fig. 
The corresponding extremal probability paper easily constructed, 
with the ordinate given and the abscissa scale given directly 

Extremal Type Distribution (Also Termed Fretchet Fisher-Tippett 
Distribution).—The Type distribution often used the prediction extreme 
wind speeds following the work Thom (33). Its application extreme wave 
prediction deriving from the use wave hindcast formulas wind speed 
distributions has subsequently been proposed [Thom (34,35), and reviewed 
Yang, Tayfun, and Fallah 


TABLE 2.—Scale Relationships for Probability Papers 


Ordinate scale Intercept 


Distribution 


(1) (3) (5) 


Type 
(Lower Bound) 
Type 
(Upper Bound) 


Fig. 2(c) shows Type density functions for different values The 
distribution corresponds In(H) possessing Type distribution, and thus 
plots straight line log-extremal paper (sometimes termed 
paper), which the ordinate scaled with extremal paper, —In [—In(P)], 
while place plotted linearly along the abscissa. 

Extremal Type III Distribution (Also Termed Weibull Distribution).—There 
are two alternative forms which have been adopted proposed for extreme 
wave prediction. These are lower-bound and upper-bound distributions, denoted 
here Types III, and which depend whether the location parameter, 
used described lower-bound (Type III,) upper-bound (Type 
the wave heights. Thus, either case represents limiting value 
the variate, beyond which occurrence possible. The Type III, 
distribution was initially developed Gumbel (14) for estimating worst possible 
droughts. Since the intention provide empirical fit the data, the 
use either Type III, justified, even though the limits the 
distribution’s shape may may not physically appealing. Figs. 2(d) and 
2(e) show typical plots linear scales these two distributions for different 


Abscissa 
Qa a 
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values For particular parameter values, the Weibull distribution reduces 
(when 1). 

will seen that the Type III distributions depend three parameters, 
rather than two with the other distributions considered, and provide added 
flexibility attempting fit data. Against this the fact that the theoretical 
variation log-extremal paper does not fall straight line. order 
obtain straight line plot becomes necessary adopt scales which depend 
one other the parameters estimated. Thus, one may plot In(H 
the abscissa, while the ordinate scale appropriate function only 
given Table Alternatively, one may plot for Type 
for Type along the ordinate, while plotted 
directly along the abscissa. Scales providing straight line plot thus depend 
making estimate before plotting the data. 

However, published applications, authors have invariably adopted chosen 
values (e.g., see Ref. (e.g., see Ref. 27), rather than best estimates 
these parameters. This simplification enables data plotted directly, 
and the Type III, distribution has been used quite successfully this way. 

the Type III, distribution, may considered correspond small 
wave height representing low level background wave activity which always 
present, lower bound the wave heights included the data sample. 
The Type distribution has not generally been adopted for extreme wave 
prediction, although its use has been advocated St. Denis (30), and Borgman 


(5). The upper height limit entailed the Type distribution makes 
particularly attractive, e.g., locations where height may limited water 
depth, fetch, other such feature, any case where finite maximum 


Data 


Various methods have been used assemble sample wave statistics 
which the extrapolation procedure applied. The individual data points 
used may defined quite different ways have been reviewed and compared 
Nolte (23). 

common method collecting wave data consists recording the free 
surface elevation over recording period (of say min) with successive 
recordings commencing the start each recording interval (of say several 
hours). The procedure for converting individual wave record into characteristic 
wave statistics well known and has been described, e.g., Tucker (36). 
For data obtained this way, one may assign one statistic each successive 
recording interval, characterizing that particular sea state, and this may chosen 
the significant wave height, expected maximum height, most probable maximum 
height, zero-crossing period, etc. This approach has been used Draper (8), 
Draper and Driver (10), Dattari (7), and others. 

common form presenting such wave data with bivariate histogram 
Each recording interval has associated significant wave 
height and (e.g., zero-crossing peak) wave period. the scatter diagram, 
all such wave heights and periods are divided into equal ranges, and the number 
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occurrences different class combinations these are tabulated. This data 
may then fitted long-term distribution wave heights. 

alternative procedure, such data may used obtain the long-term 
distribution individual wave heights, This distribution may 
given terms the long-term distribution significant wave heights, 
and the Rayleigh distribution, which describes the short-term distribu- 
tion individual wave heights, corresponding particular H,. Pedersen 
(26) expresses this distribution, 

0 
Battjes (2) proposed improvement this formula account for the correlation 
between and the characteristic wave period pertaining any recording 
interval. The modified expression 


dH, 
0 0 


which the long-term average number waves per unit time; and 
the long-term bivariate probability density and 

Nolte (23) has illustrated how unreliable design wave estimate can when 
based single year’s wave records. example provides shows 100-yr 
wave estimates differing factor about when based wave records 
two successive years. This emphasizes the importance adequate data 
base for reliable predictions made. 

The alternative procedure involving hindcasting from meteorological data has 
the considerable advantage that data deriving from relatively long time span 
may reconstructed. Such alternative commonly used for hurricane 
generated waves (e.g., see Ref. 6). statistic that generally used this 
context the most probable maximum height occurring individual storm 
greater than certain intensity [e.g., see Durning (11), Haring and Heideman 
(16), Petrauskas and Aagaard Jahns and Wheeler (18) describe the estimation 
this statistic, and also other relevant storm parameters based hindcast 
data. 

the case hurricane generated waves, Ward, Borgman, and Cardone 
(27) have described alternative approach estimating the probability distribu- 
tion large individual wave heights. This obtained terms empirical 
probability distributions hurricane tracks, affecting the closest distance 


approach the site question, hurricane intensities and the frequency 
hurricane occurrence. 


order plot the data one must first assign value P(H) each 
value the sample. this, the data are first ordered according wave 
height. The suffix then used denote its rank, with corresponding 
the largest value, and the smallest value sample containing 
wave heights. 


100 MAY 1981 
heights then given as: 


those cases which large (such when scatter diagram 
individual wave height distribution used), may calculated for fewer 
heights selected convenient intervals obtained directly from scatter diagram, 
selected the largest heights the sample. this case, the formula 
Eq. still applies, but now takes the following form: 


Number height values 


Gumbel (15) has shown that the expected probability for the mth ordered 
observation given m/(N noted the foregoing, and that this 
independent the distribution adopted. However, Eq. not unique, and 
fact has been demonstrated [e.g., see Kimball (21), and Gringorten 
introduce slight bias peculiar the distribution being estimated. unbiased 
formula depends the particular distribution being used, and also the specific 


parameters involved. more general plotting formula may written the 
form: 


which and constants determined. For the Type distribution, 
Gringorten (12) gives 0.44, 0.12; and for the Type III, distribution 
with prescribed Petrauskas and Aagaard (27) give 0.30 0.18/a, and 
0.21 0.32/a. Even so, most authors have adopted Eq. plotting 
formula for reducing the data, because its simplicity, but spite resulting 
possible departure from any underlying distribution (Ref. 12). 


Parameter 


Each the five distributions considered here actually family distribu- 
tions, whose properties depend upon the parameter values assigned them. 
Having selected one distribution likely model, remains estimate the 
parameter values which will provide the best empirical fit between the distribution 
and the data. 

The most straightforward approach plot the individual data points 
the selected probability paper and then draw straight line through these 
eye. (In the case the Type III distributions, this would apply only 
were chosen advance.) Alternatively, the best-fit line may derived 
three possible methods: (1) The method moments; (2) the method least 
squares; and (3) the method maximum likelihood. Each these three methods 
may give different estimate the parameters based the same sample. 
the following, the estimated values the parameters are distinguished 
caret over the symbol. 


Method Moments.—The method moments operates equating the first 


N+B 
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two three moments the distribution those the data. Since the lowest 
moments tend exert the strongest influence the shape distribution, 
this procedure often leads acceptable model. 

The estimated values the parameters, and are expressed terms 


TABLE 3.—Parameters Distributions Estimated Method Moments 


Distri- Estimated Parameters 
bution 


(1) 


[A - (A)’} 1/2 


FIG. 3.—Variation with Shape Parameter, for Type Distribution 


directly from the data and are defined follows: 


The log-normal, Type and Type distributions each involve only two parameters, 
and these can obtained equating the first and second moments the 
sample those the model. the case the Type distribution, the 


normal 
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relationship between and used the Type distribution sketched 
the Gamma function, and seen increasingly insensitive 
becomes large. The Type III distributions involve three parameters, 
and these, may first estimated equating the skewness, 
the sample that the model. The remaining parameters can then estimated 
equating the first and second moments the sample and the fitted distribution. 


are the second and third central moments the distribution. The sample skewness 
given 


a a Qa a 
a a 


which for the Type III, distribution; for the Type 
distribution; and the Gamma function. Thus, estimated Eq. 
which the sample skewness, Eq. applied. This result for the Type III, 
distribution was given Gumbel The relationship given Eq. plotted 
Fig. noted that the skewness becomes increasingly insensitive 

Method Least Squares.—The intention provide straight line fit 
the data when plotted the pertinent probability paper. the cases 
the lognormal, Type and Type distributions, the scales are known 
priori and only two parameters are needed each case. Thus, the least-squares 
method its most basic form directly applicable these distributions. This 
gives the slope and intercept the best-fit line terms 


FIG. 4.—Variation with Shape Parameter, for Type Ill Distributions 


ww2 LONG-TERM DISTRIBUTIONS 103 


the coordinates (x,,y,) all the data points. The corresponding estimated 
values the parameters the distribution, required, may then obtained 
from the slope and intercept the expressions given Table 

the case the Type III distributions, either the ordinate abscissa 
dependent upon one the parameters estimated, and thus without prior 
knowledge this parameter the data cannot even plotted [unless 
chosen take some trial value, carried out Battjes (2) and Petrauskas 
and Aagaard (27)]. any case, the least-squares procedure can readily 
involves iterative procedure which entails serious difficulty. 

Method Maximum Likelihood.—The method maximum likelihood attempts 
provide estimated parameters which would give the data sample the highest 
probability being observed its particular form. The available sample 
considered consist series independent observations from the same 
distribution and the probability the intersection these events the product 
their individual probabilities. The likelihood function defined as: 


which the probability density the chosen distribution; and 
the measured wave heights. The method maximum likelihood then selects 
values each parameter which maximizes [Since most the common 
density functions have exponential form, this procedure simplified 
maximizing instead In(L).] 

The method results estimated parameters which are unbiased and have 
arelatively small variance. However, the solution usually requires lengthy iterative 
manipulation, and the method consequently rarely used the present context. 
The method has been applied Thom (33) the Type distribution for 
the analysis extreme wind speeds, and its possible application the distribu- 
tions considered here described Mackenzie (22). 


Once distribution has been fitted set data one other method 
described, becomes desirable appraise the cioseness fit the data 
points the fitted distribution. The scatter the data may best described 
terms confidence limits either side the fitted line. Thus, curves 
drawn either side the best-fit line provide series confidence bands 
which indicate the confidence attached any particular data point. This 
sketched Fig. which shows, e.g., the 50% and 90% confidence bands, 
within which data are expected lie with probabilities 50% and 90%, 
respectively. (In this figure, the symbols shows example pair 
calculated height limits used derive the confidence bands.) 

procedure for constructing confidence bands was developed Gumbel 
(15), and has been described St. Denis (29) the context wave prediction. 
The method simple use, but based approximation which really 
only suitable for data near the center the ordered sample, and becomes 
invalid for the highest points the sample (which greater concern 
the present case). Furthermore, the method described applies only the Type 
distribution. 
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The more complete derivation confidence bands, applied any chosen 
distribution, described Borgman (3) and Gringorten (13). For given sample 
size rank and chosen confidence probability level, this method provides 
pair height limits for the mth statistic. These values may then plotted 
either side the mth data point, and faired line drawn through equivalent 
limits for the remaining points. This procedure indicated Fig. 
emphasized that this method cannot used for prediction (cf. Ref. 29), but 
strictly applies only the available data points. 

alternative approach has been used Petrauskas and Aagaard (27) and 
involves Monte Carlo simulation generate random sets data deriving from 
the best-fit distribution that has been obtained. The spread this simulated 
data can then used describe the confidence uncertainty attached 
any chosen value. The major advantage this approach that can used 


50% timits 


- 
- 
© 

a 


(m) 


FIG. 5.—Wave-Height Exceedence Diagram Showing Use Confidence Bands 


for any values height, including predicted values, and thus, enables confidence 
limits attached any chosen design wave height. 


Wave 


The processes selecting particular distribution, together with estimating 
the most suitable parameters, have far been described, and now consider 
methods using the best-fit line predict design wave height selected 
correspond prescribed return period. The return period (or recurrence 
interval), the average time interval between successive events the 
design wave being equalled exceeded. The wave height selected for this 
purpose, would initially defined the same way the recorded wave heights 
(e.g., the significant wave height over recording interval, the most probable 
maximum height individual storm). The return period, may expressed 
dimensionless form the expected waiting time, which 
the recording interval associated with each data point. (Thus, the case 


.00001 
os 
i. 
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where each data point corresponds individual storm, would the average 
time interval between such storms.) The expected waiting time directly related 
the probability, P(H), 


From this formula, return periods may scaled along the ordinate, indicated 
Fig. Thus, prescribed return period has associated value 
and the corresponding wave height may then determined from the (extrapolated) 
best-fit line that has been plotted. This design value will defined the 
same way were the individual data points (e.g., the significant height over 
recording interval), and should equalled exceeded the average once 
over the duration any return period. The probability distribution individual 
wave heights during the corresponding extreme sea state, other characteristic 
heights that sea state, such the most probable maximum height the 
height corresponding chosen probability value, may then obtained 
terms the design height obtained the use statistical formulas applicable 
short-term (stationary) variations. 

easily shown that the return period for the common ordinate scale, 
used the Types II, and distributions may 
approximated 


Nolte (23) has indicated that many cases correction the return period 
necessary, because its intended meaning distorted due the occurrence 
more than one data point single storm; i.e., the occurrence large 
valued data points will tend group together, and the average duration between 
successive storms containing design wave, which really the intended 
requirement, will thus larger than the true return period. The corresponding 
correction factor’’ relatively important when individual wave 
model used (since more data points will tend occur groups), and 
unnecessary when storm model used 


many cases, design carried out with wave heights corresponding 
some prescribed return period, e.g., 100 yr. However, appropriate 
the context offshore design consider also the encounter probability, 
This the probability that the design wave equalled exceeded during 
prescribed period, say the design life-time structure; and indeed 
may preferable instead select correspond prescribed values 
and The relationship between these quantities (e.g., see Ref. is: 


r L/r 


which generally when T/r >>1), suitable approximation 
which independent the recording interval is: 


Thus, e.g., the return period giving rise encounter probability 0.1 for 
design life-time yr, 475 yr; and again, the probabilities 50-yr 
design wave being exceeded within yr, yr, and 100 are 0.181, 0.632, 
and 0.865, respectively. 

the design life-time the structure prescribed, the encounter probability, 


may scaled along the ordinate alongside, place of, sketched 
Fig. 


Extreme Wave 


There are various approaches estimating the extreme wave period associated 
with the design wave height which has been obtained. The first repeat 
the entire procedure using wave periods instead wave heights statistics. 
Thus, when using scatter diagram, marginal frequencies are obtained summing 
the number wave occurrences each period range. using the same return 
period the heights, one may obtain predicted value the 50-yr zero-crossing 
period, for hypothetical wave record with the same recording interval. 
Draper (8) has suggested that this value may used with the predicted 
height, but this suggestion based the fact that there noticeable correlation 
between the two variables the scatter diagram. 

alternative approach involves using the predicted wave height set 
lower limit the wave period (due wave breaking), and then using series 
values the wave period above this lower limit find the worst effect 
the structure. assuming Pierson-Moskowitz spectrum, Battjes (2) has 
shown that for deep water and intermediate depths, characteristic wave steepness 


which the gravitational constant [9.801 m/s’ (32.2 Thus, 
lower limit period for given significant height, may set as: 


Research has recently been directed Houmb and Overvik (17), and 
Ochi (24)] towards fitting bivariate probability distribution long-term data 
both wave heights and periods. This predicts different combinations height 
and period corresponding known probabilities occurrence and should account 
for limits height-period combinations due wave breaking. Such distribution 
can conveniently presented contours known probabilities scatter 
diagram. Ochi (24) has employed bivariate log-normal distribution this way. 


APPLICATION 


variety examples the foregoing procedures have been provided 
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the various references cited herein and will not repeated here. These have 
invariably adopted one distribution another fit the wave data. 
clear, though, that estimated design wave height from given set 
data depends some extent the particular distribution selected and the 
method curve fitting applied. Moreover, many cases (e.g., see Ref. 7), 
the tail the data, corresponding the most extreme measurements made, 
shows the most noticeable departure from curve fitted all the data, even 
though this extreme range precisely that over which the best possible fit 
required. Thus, may preferable apply, e.g., least-squares fit 
the tail the data rather than all it. 

This approach has been explored MacKenzie (22). computer program 
has been prepared the writers which wave height data analyzed 
predict design values. their approach, any one probability distribution 
not chosen priori, but rather the five distributions described herein are all 
fitted the data using least-squares fit—a three-parameter fit used 
the case the Type III distributions. Confidence bands are then constructed 
about the fitted curves and the basis these the most appropriate distribution 
selected. This then used derive design wave heights the usual way. 
This procedure has been illustrated MacKenzie (22) using data gathered 
Tinier Point, New Brunswick. 


Summary 


The various steps involved the estimation extreme design waves have 
been reviewed. These include methods data collection, the use plotting 
formula, the selection suitable distribution and its parameters, the plotting 
confidence bands about the best-fit line, and the selection design wave 
corresponding prescribed return period encounter probability. 

The probability distributions that have been used proposed include the 
log-normal and Extremal Types II, and III distributions (Gumbel, Fretchet, 
and Weibull distributions, respectively). These distributions are described and 
their properties and other results required for fitting data are assembled 
tabular form for ready reference. 

Additional research required enable bivariate probability distributions 
fitted long-term wave height and period data, taking into account the 
limits imposed wave breaking. 
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The following symbols are used this paper: 


slope; 

intercept; 

encounter probability; 

gravitational constant; 

wave height; 

wave height mth data point; 

significant wave height; 

long-term cumulative probability individual wave heights; 
design life-time; 

rank; 

sample size; 

cumulative probability; 

probability density; 

cumulative probability Rayleigh distribution; 
recording interval; 

wave period; 

long-term average number waves per unit time; 
return period; 

zero-crossing period; 

abscissa scale; 

ordinate scale; 

shape parameter; 

skewness; 

gamma function; 

location parameter; and 

scale parameter. 
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ASCE. This paper part the Journal the Waterway, Port, Coastal and Ocean 
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No. WW2, May, 1981. 
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TECHNICAL NOTES 


CABLE USING ORTHOGONAL 
COLLOCATION 


Herman and Ernest McReynolds” 


The capability analyze and design cable systems during the deployment 
and retrieval phase must encompass the dynamic effects paying out and 
reeling the cable system. Past cable payout/reel work has been prepon- 
derantly treated quasidynamic manner (4). With this quasidynamic approach, 
dynamic analysis was performed fixed length, for several lengths 
cable. More recently, spatially-discrete methods, such the lumped parameter 
and finite element methods, have been employed manner where the length 
the cable system changed each time integration step (11). These adjustments 
present rather adrupt change configuration and result, numerical 
oscillations can introduced. Alternatively, cable systems undergoing changing 
length can treated continuous fashion where full dynamic effects due 
payout reel are considered (5). 

Consider spherical load mass, and radius, supported extensible 
permeable cable; S(t) the unstretched length cable between the load and 
the top support time, Under the assumption that only vertical motion 
present, the equations motion for this one-dimensional system are: 


S(t), 


s= S(t), v [S(t),¢] 


which elastic modulus (assumed linear); cable cross section 


Prof. Mech. Engrg., Portland State Univ., P.O. Box 751, Portland, Oreg. 
97207. 

Assoc., Portland State Univ., P.O. Box 751, Portland, Oreg. 97207. 

Note.—Discussion open until October 1981. extend the closing date one month, 
written request must filed with the Manager Technical and Professional Publications, 
ASCE. Manuscript was submitted for review for possible publication June 1980. 
This paper part Journal the Waterway, Port, Coastal and Ocean Division, 
Proceedings the American Society Civil Engineers, Vol. 107, No. WW2, 
May, 1981. ISSN 0148-9895 /$01.00. 
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area; cable mass per unit unstretched length cable; mass 
displaced fluid per unit length cable; added mass load the fluid 
(the sphere); mass displaced fluid the load (the sphere); C,, 
drag coefficient load; drag coefficient tangential cable; fluid 
density (per unit volume); diameter cable; and deterministic 
ship motion. 

The independent variables are the time, and the material coordinate, 
which the distance along the unstretched length from the load the cable 
mass element. The dependent variables, and the tension and vertical 
velocity (positive up), respectively. Eqs. and are the usual equations for 
longitudinal oscillations elastic cable. Boundary conditions are given 
and 

The boundary condition Eq. very awkward handle since applied 
variable material coordinate, i.e., S(t) which S(t) the unstretched 
length cable the system. resolve this difficulty, material coordinate 
transformation will made that the new not time 
dependent. This accomplished the following change variables: 

The transformed equations are: 

S(t) 


m—m 


The method weighted residuals, particular that orthogonal collocation, 
provides approach approximate the solution the preceding nonlinear 
system partial differential equations. 

positive integrable weight function, chosen integrating weight 
obtain Gaussian quadrature for any integrable function, for 
The quadrature formula the form 


k=l 


which the nonzero coefficients; and the values lying the interval 


denotes the number terms the preceding summation, and then 
varied over the positive integers, set polynomials obtained: 


Re +o——- 
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For the selected weight function, these polynomials will orthogonal over 
the interval integration (3). 

Trial solutions Eqs. 5-8 are taken linear combination polynomials 
from this set having degree less than equal for some integer 
Following the method Galerkin, i.e., taking the interior equations that 
system with these expansions substituted interior residuals (1): 

S(t) 


S(t) 


The inner product these residuals are taken with the first polynomials 
over (0,1) and set zero. Thus 


the preceding integrals yields: 


For each the aforementioned equations determine homogeneous system 
matrix given by: 


A= (ay), ay = P(o,) 


This matrix composed linearly independent functions evaluated 
distinct points which will make nonsingular. Therefore, the preceding system 
must have the solution 


j=1,2 


the aforementioned set equations are added the boundary conditions 
The set points (0, are the so-called orthogonal collocation 
points (9). 

Derivatives with respect the normalized coordinate appearing the residuals 
are eliminated variables the consideration the polynomial trial solutions 
(2). For example: 


(16) 


116 MAY 1981 


Differentiating and evaluating this quantity the collocation points yields 


n+l av n+l 


functions evaluated distinct points which will allow the aforementioned 
obtained. These equations were solved using the Hamming predictor-corrector 
method with fourth order Runge-Kutta method startup procedure (8). 

Available experimental data were then compared analytic results order 
gain confidence applying the developed computer program (10). The model 
anchor was sphere. Two different types cable were used; they were solid 
rubber and braided nylon. The cable and anchor were reeled paid out 
the top tank means winching mechanism. The winching mechanism 
also provided measure the payout/reel velocity and acceleration. The 
simulation ship motion was accomplished oscillator which caused the 
winching mechanism and winching platform oscillate the vertical direction. 
Cable tension was measured the top. 

Preliminary comparison between experimental and analytic results revealed 
the importance modeling the payout/reel acceleration. Spline functions 
were employed continuously approximate the experimentally-derived data, 
and these conditions were input into the collocation program. For both experi- 
mental and analytic results, the response the modeled system, terms 
tension this winch end, exhibited transient oscillation due the initiation 
payout/reel for both experimental and analytic results. The tension then 
settled steady-state configuration encompassing effects due change 
system length and hydrodynamic drag forces the cable. The salient features 
the analytic results agreed reasonably well with the experimental data; however 
amplitude comparisons revealed significant differences, that the computer 
output tended overactive. 

Two remedies were incorporated into the computer program. Dynamic uniaxial 
tension tests revealed measurable material damping for the case the rubber 
cable. Material damping was incorporated the analytic model for rubber (6). 
The other problem was specification system drag the case braided 
nylon. The system was subjected velocity conditions which were successively 
transient, periodic, and constant nature. Based sophistication involved 
with conducting previous drag studies (7), attempt could made conclu- 
sively investigate drag for this particular system. result, load drag was 
somewhat arbitrarily increased from for the nylon case. Steady-state 
response from the experimental case revealed constant component for increased 
drag, irrespective parametric variation EA, masses, cable drag. 

Since the same payout/reel program was utilized for parametric analysis, 
numerical difficulties cannot ruled out contribution. Lower order numerical 
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integration was employed due word-size limitations and the procedure for 
approximating derivatives. 


The method weighted residual, orthogonal collocation technique was used 
solve the continuous equations motion for cable payout and reel in. The 
method was compared experimental data for two types cable under payout 
and reel with and without ship motion excitation. Results were within about 
35% for the case rubber cable using constant coefficient drag. For 
the case nylon cable, the collocation method showed similar results relative 
experiments, but some differences phase and frequency were encountered. 
Changes coefficient drag lessened these differences. demonstrated 
both analytic and experimental results significant dynamic tension were 
encountered due solely initiating and bringing the cable system 
abrupt stop. 

suggested previous sections, the thrust this investigation was develop- 
ment and exercise numerical solution technique for accepted cable 
model (6). Future work will include alternative numerical integration procedures 
and investigation NWR solution cable-fluid models which may not in- 
corporate velocity-squared drag. addition, two-dimensional solution the 
current form equations motion are under development and will compared 
experimental results. 
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Closure Philip Vitale* 


The writer thanks the three discussors for the significant and interesting 
comments they have made. Eq. 31, derived Jonsson, concise relation 
which will prove most useful. Lofquist’s use Eq. relate Eqs. and 
provides valuable link between different sections the report. Both 
these analyses suggest that the factors the rippled bed trends 
may due simply the constant period the data. Nielsen’s Eq. 
extends the original results and allows for more direct comparison with Jonsson’s 
values The relationship between and are discussed Refs. and 
15. 

point brought all three discussors concerns the small contribution 
the ripples the total resistance, predicted 25. Lofquist notes 
that the sand grain component the roughness (see Eq. 21) includes some 
effect the ripples due the alteration the flow patterns and skin friction 
distribution; therefore, Eq. does not take the into account 
completely. Nielsen mentions two references (21,22) which demonstrate the 
importance sand grain roughness. Jonsson suggest that the 30% contribution 
the ripples too small. ‘his point, the writer refers his interpretation 
the data which seemed best explained using Eq. for the bed roughness. 
This can seen the series graphs Figs. 10-13. Fig. divides the 
data into three groups: the turbulent data, defined Fig. and represented 
the symbol the Area data, defined Fig. and represented the 
symbol and the balance the data which represented the symbol 
Fig. 11, which uses Eq. for the abscissa, the rippled and flatbed 
data group together allow single straight line drawn through the 
data. the ripple effect increased the abscissas Figs. and 13, 
the rippled data move away from the flat-bed data making necessary draw 
one line for each group. 

The writer recognizes the problem with comparing Eqs. and past 
determinations 2.5 for flat beds and for rippled beds. 
fact, was this existence two different relations which led the writer 
the analysis the data Carstens, Neilsen, and Altinbilek (4) attempt 
find roughness factor which would allow for smooth variation the 
roughness the bed changes, with increasing flow, from rippled flat-bed 
sheet flow. (Lofquist points out that the flat-bed condition which Eq. 
applies sheet flow.) Eqs. and allow for this smooth transition. 

August, 1979, Philip Vitale (Proc. Paper 14754). 


Hydr. Engr., Dept. Army, U.S. Army Coastal Engrg. Research Center, Corps 
Engrs., Kingman Building, Fort Belvoir, Va. 22060. 
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The issues the preceding two paragraphs inspire additional tests which 
would specifically look into the relative roles ripple and sand grain roughnesses. 
These would include tests similar fixed-bed smooth, fixed-bed rough, and 


10.—Data Separated into Three Groups Defined Text 


(2 +0.01 


FIG. Abscissa; Rippled and Data Group Allow Single Straight 
Line 


movable-bed ripples under identical flow conditions. Also, flow over movable-bed 
ripples can gradually increased until sheet flow produced while the roughness 
factor continually determined track its evolution. 

regards the turbulent zone Fig. the writer was purposefully conservative 
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placing insure that the points used Fig. were indeed turbulent. 
This demonstrated the tight grouping the data points Fig. 11. 
possible make good case, Jonsson does, for shifting the line 


FIG. 12.—Ripple Effect Increased 


FIG. 13.—Greater Increase Ripple Effect 


the left. This would then include the turbulent zone number the 
data points Fig. 11, some which fall close the group. 
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The writer would like thank Young and Sobey, and also Venezian for 
their helpful comments. Both discussions relate practical alternatives iterative 
solution the wave dispersion equation. 

Young and Sobey draw attention the tabular procedure, and compare its 
computational efficiency with the writer’s approximate solutions, with very 
favorable results. Some the gain probably arises from exploiting the deep-water 
simplification for but some advantage clearly apparent when sufficient 
computer storage available. 

Venezian raises the important question balance between range validity 
and accuracy, and presents two shallow water Padé approximations, Eqs. 
and 13, which are highly efficient the range These expressions 
are not the only approximants that could used this range, nor indeed 
was the writer’s Eq. the only approximant that could used 

Both discussions emphasize the desirability subdividing the range 
and may that more than two subintervals would further improve 
computational efficiency. There remains the problem balancing computer 


time, computer storage, and simplicity programming, for given computational 
accuracy. 


Errata.—The following corrections should made the original paper: 


Page 457, line Should read the relation instead 


Page 458, paragraph line Should read Ref. instead 
Ref. 


Page 458, paragraph line Should read term instead 


“November, 1979, John Hunt (Proc. Paper 14944). 


Applied Mathematics, Univ. Reading, White Knights, Reading RG6 2AX, 
England. 
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IRREGULARITIES SOLUTIONS NONLINEAR WAVE 
VERTICAL 


the original paper the author pointed out Isaacson’s incorrect conclusion 
(24) that perturbation solution for the diffraction problem may not formally 
extended beyond the linearized approximation due the irregularity higher- 
order solutions along the contour intersection the vertical pile and the 
free surface, and suggested procedure compute the solution the contour 
intersection. 

Although the present contributor agrees with Miloh’s view the incorrectness 
Isaacson’s conclusions, some further consideration needs paid before 
the adoption the solution proposed the original paper. The aim this 
contribution therefore complement the views Miloh shedding more 
light into the origin and the nature these irregularities which may then lead 
solution procedure. 

Since the apparent ease using polar cylindrical coordinate seems hinder 
the source this irregularity, shall use different notation. Then the boundary 
conditions creating the irregular behavior may written as: 


and scattered potentials; subscripts and denote the derivation 
with respect time and normal, respectively; and the sets consisting 
the points the free surface and body surface, respectively. 

the points intersection form new set the first question answered 
the choice between the following options: (1) S,, (2) 

For the first second option there problem and thus irregularity 
the solution, since the usual perturbation solution without any correction 
satisfies the needs these problems. If, however, the third option the real 
situation, then the problem not simply the satisfaction Eqs. and Rewriting 


(k) 


can detect that: 


1980, Touvia Miloh (Proc. Paper 15385). 


Appointment, The British Ship Research Association, Wallsend Research Station, 
Wallsend, Tyne and Wear, NE28 6UY, England. 
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create irregularity, solution which will normally require the extension 
two-dimensional complex variable techniques three-dimensional. 

For points belonging S,, the direction normal undefined. the 
problem treated the author this manifests itself radial vertical derivatives 
depending the direction approach. study the boundary discontinuity 
effect for Neumann problem, Craggs, Mangler and Zamir (20) indicated that 
the behavior the singularity strongly depends the vertex angle. 


Since Rankine source distribution both the body surface and the free 
surface also leads solution, done Gadd (23) and Dawson (21) 
the ship-wave resistance theory, the simultaneous satisfaction Eqs. and 
first raises the question the existence solution since these 
conditions may overdetermine the system equations. 

The remedy the second singularity arising from the nonuniqueness the 
boundary normal may easier treat. Recently Eggers (22) gave rigorous 
treatment the singularities arising from the boundary discontinuities which 
may extended the present case. 

From the foregoing discussions clear that neither the problem nor the 
origin the irregularity well-defined unless one the options the orientation 
the set chosen. Furthermore, there irregularity this would also 
exist the first-order solution since both the jump the coefficients and 
the nonuniqueness normal also exist the first-order problem. 


20. Craggs, W., Mangler, and Zamir, M., Remarks the Behaviour 
Surface Source Distribution Near the Edge Aeronautical Quarterly, 
Vol. 24, Part 1973, pp. 25-33. 

21. Dawson, W., Practical Computer Method for Solving Ship-Wave 
Proceedings, Second International Conference Numerical Ship Hydrodynamics, 
1977, pp. 30-38. 

22. Eggers, K., Irregularities the Wave-Flow Due Source Panels and How 
Compensate them Adding Kelvin-Source Line Elements,’’ Contribution the 
Oct. 10-12, 1980, Continued Workshop Wave Resistance Calculations, held 
Izu, Japan. 

Gadd, E., Method Computing the Flow and Surface Wave Pattern Around 
Full Transactions the Royal Institute Naval Archives, Vol. 118, 1976, 
pp. 207-219. 

Isaacson, St. Q., Wave Forces Large Offshore 
Journal the Waterways, Harbors, and Coastal Engineering Division, ASCE, Vol. 
103, No. 1977, pp. 166-170. 
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Errata 


The following reference should added the original paper: 


Thompson, F., Spectra Shallow U.S. Coastal 
Technical Paper No. 80-2, U.S. Army Coastal Engineering Research Center, 
Fort Belvoir, Va., Feb., 1980. 


INLINE Forces VERTICAL CYLINDER Waves” 
Discussion George Christodoulou,’ ASCE 


The author has presented noteworthy experimental study concerning the 
inline forces acting fixed smooth vertical cylinders wave tank. Detailed 
force measurements and application fifth-order stream function wave theory 
resulted defining relationship between the (local) Morison equation coeffi- 
cients C,,, and the Keulegan-Carpenter number, (Fig. 3). The potential 
value the study improving offshore structure design obvious. 

However, there severe limitation the present results, associated with 
the Reynolds number parameter. Besides its small values, distinction 
could practically made, because the small range variation, i.e., 
10*-3 The author suggests that the scatter the experimental data 
about the mean curves Fig. partly due the change Previous 
studies (3,5,15) indicate strong dependence the coefficients, especially 
Therefore may further suggested that the Reynolds number change 
responsible, some extent, for the trend the mean curves well, i.e., 
the gradual decrease and increase C,, for 10. This would 
consequence positive correlation between and KC. Such remark 
supported Fig. where observed that the curves C,,, obtained 
the author are generally closer the broken lines for 


1980, Edward Thompson (Proc. Paper 15385). 
1980, Subrata Chakrabarti (Proc. Paper 15403). 
Applied Hydr. Lab., National Technical Univ. Athens, Athens 624, Greece. 
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can further seen Fig. that there appreciable difference between 
the values the coefficients C,,, obtained for wave-flow and those for 
two-dimensional harmonic flow. The major discrepancies occur the high end 
the range, where the sparsity data makes the wave-flow results less 
reliable, the author has mentioned. The agreement remarkable view 
the empirical nature these coefficients, their dependence the wave 
theory describing the flow field, and some unrealistic assumptions commonly 
involved their evaluation, e.g., constancy over the wave cycle. Certainly, 
wave flow experiments are more representative field conditions than harmonic 
flow ones, however the limited present results lead the tentative conclusion 
that significant differences the coefficients should expected for the 
two cases. This could clarified only further similar experiments, which 
would cover wider range and, most importantly, extend higher ranges 
encountered the field. would very useful this procedure could establish 
more reliable estimates the Morison coefficients, the lack which necessarily 
leads oversimplifications design practice (15). 

Nevertheless, the more accurate evaluation the Morison coefficients may 
not the only way proceed future research. Given the drawbacks associated 
with their empirical nature and the immense growth computing capabilities 
since the Morison formula was proposed, more generalized correlation 
wave characteristics with resulting forces seems now both feasible and more 
sound. For example, the approach proposed Jen (14) could lead such 
generalized relationships. 

Concerning the comparison between measured and calculated total forces 
the cylinder, shown Figs. should kept mind that there 
considerable between the forces measured and those calculated. 
Indeed, the coefficients used the theoretical force calculation had been obtained 
local force measurements the same cylinders, the same experimental 
arrangement and using, backwards, the same wave theory, Therefore, one 
would expect good agreement between essentially the local forces and the 
integrated forces the same members. may noticed that, although still 
satisfactory considering the complexity the problem, the agreement relatively 
poorer for the shallow cylinder (Fig. 6), which the one with the smallest 
input the evaluation the mean C,,, curves. better establish the 
validity the proposed functions C,,, versus would seem essential 
compare them additional data derived from independent experiments. 


14. Jen, Y., Force Analysis: Alternate Journal the Waterway, 
Port, Coastal and Ocean Division, ASCE, Vol. 106, No. Proc. Paper 15166, 
Feb., 1980, pp. 117-121. 

15. U.S. Army Corps Engineers, Shore Protection Manual, Coastal Engineering Research 
Center, Vol. II, 1975. 
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Discussion Glendon Stevens, Jr.* and Claude 
Members, ASCE 


After reviewing the original paper the writers would like state that the 
authors have done commendable job constructing broad array data, 
alluding the impact gravel mining river stability. However, certain 
facets the data and the conclusions drawn therefrom inspire comment. 

the outset, the title this article appears misleading one carefully 
examines the content. becomes readily apparent perusing the body 
the article that the authors are discussing descriptor, armoring 
is, the very best, vaporous term. Simply put, refers the natural 
phenomenon prevention erosion scour the banks bed alluvial 
river. The term derived from supposition that the fines, sand and silt 
the bed, are placed suspension scour, leaving materials which are too 
large carried the flow. That supposition fails take into account 
some factors which are difficult ignore. First all, the armoring concept 
presented might plausible if: (1) There was increase flow sufficient 
increase the carrying capacity; and (2) there was inflow fines from 
upstream. Such conditions might exist below large reservoir. However, 
postulate that gravel the size normally mined for commercial purposes can 
prevent scour river subjected wide range flow the Mississippi, 
dance the fringes fantasy. After all, concrete mattresses are being 
destroyed constantly. 

The authors have unrestrainedly discussed the detrimental effects commercial 
gravel mining, ignoring completely the possibility that the dredged areas could 
serve storage areas within which the sediment naturally carried the river 
could trapped. Sediments trapped this manner are unlikely hazard 
navigation since most commercial gravel mining takes place well aside from 
the navigation channel. the process spoken operates proclaimed, then 
should alleviate small degree the need for channel dredging, regard 
maintenance open-channels for navigation. 

addition substantive questions reference this article, this reviewer 
question the location the photographs shown Figs. 
and These pictures seem reveal that the lower Mississippi River was 
dry the time they were taken. The history the Mississippi River records 
such condition. The writer believes the gravel found the lower Mississippi 

1980, Peter Lagasse, Brien Winkley, and Daryl Simmons (Proc. 
Paper 15643). 


Prof. Civ. Engrg., Univ. Missouri-Rolla, Inst. River Studies, 302 
Engrg. Research Labs., Rolla, Mo. 65401. 


St. Louis Dist. Corps Engrs., 210 12th Street, St. Louis, Mo. 
63101. 
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River, because years transport and deposition, would smoother and 
more rounded than that depicted. (Reproduction might have caused distortion.) 

should noted critique this article that the amount material 
commercially dredged (sand and gravel) from the Mississippi River during any 
given year less than the material naturally transported the river. 
Therefore, the authors are correct when they state that some all the 
changes examined the paper might have occurred without the removal 
gravel from the river bed. 

addition, the writer compelled point out that when data are collected 
and manipulated the manner they were plotting Figs. and they 
lose their identity. Data collected specific point analyzed this fashion 
may acceptable. However, data collected over 200-mile reach under varying 
flow conditions and analyzed they were questionable. 

Further, Fig. indicates that relatively few years the bed material 
the river will achieve uniform gradation. This proposition which the 
writer has yet find agreement and one which exceedingly difficult accept. 
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